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'  .■  SECTION  1 

INTRODUCTION  1  1 

1  I  i 

i 

Since  the  conception  in  1959  of  using  a  large  underground 
cavity  to  decouple  (and  thus  disguise  or  conceal),  a  nuclear  test 
(Reference  1) ,  several  experimental  and  theoretical  programs  have 
been  directed  toward  improving  and  extending  this  concept.  The 
results  of  several  theoretical  studies  of  large-cavity  decou¬ 
pling  phenomena  (Reference  2)  and  the  data  obtained  from  both 
large  and  small  field  experiments  (References  3  and  4)  have  all 
served  to  confirm  the  validity  cf  the  initial  decoupling  theory, 
The  results  of.  some  small-sc^le  laboratory  experiments  conducted  ; 
on  cavities  in  plaster  of  paris  (Reference  5)  indicate  however 
that  the  theory  of  large-cavity  decoupling  may  pot  be  entirely 
valid.  1  i  ! 

i 

Although  the  basic  thedry  for  large-cavity  decoupling  is 
now  well  established,  there  are  several  aspects  of  the  problem 
that  still  need  clarification.  In  particular,  the  full  effects 
of  overburden  pressure,  non-linear  material  properties,  and  cavity 
shape  upon  decoupling t phenomena  are  not  now  well  known.  Before  a 
complete  evaluation  of  the  potential  applications  and  limitations 
of  large-cavity  decoupling  may  be  made,  the  effects  of  these  ( 
parameters  must  be  better  understood. 

*  ‘  I 

Since  June  of  1968,  Physics  International  Company  has  been 
conducting  both  computational  and  experimental  studies  of  de-  , 
coupling  phenomena.  These  efforts  at  Physics  International, 


.1 


supported  by  the  Advanced  Research  Projects  Agency  through  the 
Defense  Atomic  Support  Agency,  were  initially  directed  toward 
using  finite-difference  calculations  to  model  various  shots 
within  the  COWBOY  scries,  the  SALMON  and  STERLING  nuclear  events, 
and  various  tescs  to  be  conducted  in  the  MIRACLE  PLAY  series. 

In  support  of  these  calculational  efforts,  a  study  of  the  feasi¬ 
bility  of  conducting  small-scale  experiments  was  initiated  in 
1969.  The  results  and  conclusions  of  these  calculations  and 
initial  experimental  studies  are  discussed  in  detail  in  Refer¬ 
ence  6  . 

Beginning  in  March  of  1970,  Physics  International  under¬ 
took  a  program  for  the  Advanced  Research  Piojects  Agency  (under 
Contract  No.  DASA-01-70-C-0094)  to  conduct  several  small-scale 
experiments  and  to  perform  further  calculations  i,i  support  of 
the  MIRACLE  PLAY  series  and  in  decoupling  concepts  in  general. 

This  report  discusses  the  results  of  the  small-scale  experi¬ 
ments  carried  out  on  this  program.  Despite  several  initial 
experimental  difficulties,  a  reliable  and  reproducible  means  of 
conducting  small-scale  experiments  was  developed.  Nine  small- 
scale  tests  were  successfully  conducted  using  a  variety  of  salt 
types,  initial  detonable  gas  mixtures  (CH^  and  C^) ,  and  equili¬ 
brium  cavity  pressures  after  detonation.  All  nine  shots  were 
fired  with  a  confining  pressure  of  2,500  psi  to  simulate  over¬ 
burden  effects.  Dynamic  strain  gages,  used  to  measure  radial 
displacement  as  a  function  of  time,  gave  consistent  high-quality 
data  on  decoupling  behavior.  The  experimental  results  confirmed 
the  theory  of  large  diameter  cavity  decoupling  and  demonstrated 
the  importance  of  rock  strength  independent  of  overburden  pres¬ 
sure  . 
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The  calculational  technique  described  in  this  report  is 
applicable  to  one  and  two-dimensional  calculations.  Fifteen  one¬ 
dimensional  calculations  with  detonabls  gas  (CH4  and  02)  and 
nuclear  sources  were  performed.  Yield  models,  equation  of  state 
models,  overburden  and  cavity  equilibrium  pressures  were  varied, 
zoning  schemes  and  stability  criteria  were  examined  and  response 
and  Fourier  spectral  analyses  were  carried  out  for  each  case.  A 
two-dimensional  calculation  with  a  pill  box  cavity  was  partially 
completed,  but  was  terminated  early  because  of  lack  of  funds. 

All  calculations  made  use  of  the  LEEK  boundary  condition  which  is 
a  linear  elastic  region  coupled  to  the  full  elastic-plastic  des 
cription  of  the  earth  media,  to  extend  the  calculational  grid  to 
large  distances  from  the  cavity  without  loss  of  resolution.  This 
calculational  technique  in  conjunction  with  advanced  computer 
facilities  provides  a  very  useful  and  potentially  inexpensive 
method  for  studying  decoupling,  earthquake  and  other  ground 
motion  phenomena,  at  near  teleseismic  distances  (several  kilo¬ 
meters)  while  retaining  a  detailed  description  of  the  source 
mechanism  and  the  non-linear  earth  properties  near  the  cavity  can 
affect  the  character  of  the  teleseismic  signals. 

Section  2  of  this  report  gives  the  background  of  decoupling 
work  pertinent  to  this  investigation.  The  method  results  and 
analysis  of  the  small  scale  decoupling  experiments  are  presented 
in  Section  3.  The  method  results  and  analysis  of  the  calcula¬ 
tional  phase  of  the  program  are  given  in  Section  4.  Finally,  the 
conclusions  of  this  work  are  detailed  in  Section  5. 
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SECTION  2 
BACKGROUND 


The  theory  for  large-cavity  decoupling  as  initially  pre¬ 
sented  by  Latter,  e t  al .  ,  in  1959  (Reference  1)  ,  was  restricted 
in  its  application  to  purely  elastic  media  and  could  not  be 
applied  to  media  in  which  plastic  deformation  occurred.  Except 
for  very  low  stress  level  loadings,  all  earth  materials  in  which 
a  decoupled  nuclear  test  might  be  conducted  will  invariably 
undergo  some  plastic  or  brittle  deformation  of  a  permanent  nature. 
Consequently,  it  would  not  be  possible  to  predict  the  decoupling 
factors  for  cavities  in  earth  media  subject  to  such  nonelastic 
behavior.  Some  of  the  effects  and  implications  of  brittle  and 
plastic  rock  deformation  were  considered  by  Latter,  ejt  ail .  ,  in 
1961  (Reference  7)  .  They  concluded  that  for  decoupled  shots  for 
a  low  overburden  pressure,  such  as  existed  for  the  shallow  COWBOY 
experiments,  tensile  radial  cracking  would  be  the  failure  mechan¬ 
ism  which  would  most  likely  reduce  the  decoupling  effectiveness 
of  the  cavity.  On  the  other  hand,  it  was  concluded  that,  in  a 
medium  such  as  salt  with  overburden  pressures  resulting  from 
depths  of  burial  on  the  order  of  one  kilometer,  plastic  flow 
would  be  the  failure  mechanism  reducing  decoupling  effectiveness. 
They  noted, as  might  be  expected , that  plastic  yielding  occurring 
during  construction  of  a  deep  cavity  would  not  affect  the  de- 
coupxing  effectiveness  of  that  cavity  so  long  as  the  loading  by 
an  explosion  in  the  cavity  was  purely  elastic.  It  was  also 
concluded  thac  a  material  displaying  extreme  work -hardening 
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effects,  such  as  salt  dees  under  appropriate  conditions,  would 
tend  to  provide  better  decoupling  behavior  than  a  non-work-harden- 
mg  material.  The  purely  analytic  nature  of  the  results  pre¬ 
sented  by  these  investigators  makes  it  practically  impossible  to 
derive  the  decoupling  effectiveness  of  materials  displaying  com¬ 
plex  yielding  or  compaction  phenomena.  Generally,  the  implica¬ 
tions  of  such  phenomena  can  only  be  fully  evaluated  with  the  use 
of  finite-difference  calculations  or  appropriate  experiments. 

The  first  experimental  work  to  study  and  evaluate  the  theory 
of  large  cavity  decoupling  was  carried  out  on  project  COWBOY 
(Reference  3)  .  The  17  high -explosive  test  shots  conducted  in 
the  COWBOY  series  were  designed  to  yield  fully  tamped,  fully 
decoupled,  and  partially  decoupled  conditions.  Decoupling  fac¬ 
tors  in  general  agreement  with  those  predicted  from  the  theories 
of  Latter,  et  al.  (References  1,  7,  and  8)  were  obtained.  As 
has  already  been  noted,  the  biggest  limitation  of  these  experi¬ 
ments  was  the  shallow  depth  of  burial  and  the  concomitant  low 
overburden  pressure  to  which  the  cavities  were  subjected.  Because 
the  results  were  probably  influenced  by  the  low  tensile  strength 
of  salt,  these  experiments  do  not  provide  good  information  on 
either  the  effects  of  compressive  rock  strength  or  varying  over¬ 
burden  pressures. 

In  an  effort  to  obtain  data  on  the  decoupling  theory  where 
larger  overburden  p  essures  would  be  involved  and  where  loading 
could  be  accomplished  by  a  nuclear  detonation,  the  SALMON  and 
STERLING  experiments  were  conducted  in  the  Tatum  Salt  Dome, 
Mississippi.  The  SALMON  experiment  involved  the  detonation  of 
a  fully  tamped  5.3-kiloton  nuclear  device  at  a  depth  of  2,600  feet 
(Reference  9)  .  The  cavity  formed  by  the  SALMON  event  was  very 
closely  spherical  in  shape  with  a  radius  of  17  meters  (55  feet) . 
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The  SALMON  cavity  was  used  as  the  decoupling  cavity  for  the 
STERLxNG  event,  which  involved  detonation  of  a  0 . 380-kiloton 
nuclear  device.  The  greatest  limitation  in  deriving  good 
decoupling  data  from  the  SALMON  and  STERLING  events  resulted  from 
scaling  down  the  field  measurements  of  the  high-yield  SALMON 
detonation  to  be  representative  of  the  ground  motions  which  would 
be  observed  for  a  fully  tamped  lower-yield  device  comparable  to 
that  used  in  the  STERLING  event. 

A  third  program  designed  to  obtain  field  data  on  decoupling 
phenomena  is  represented  by  the  MIRACLE  PLAY  series.  The  princi¬ 
pal  objective  of  the  MIRACLE  PLAY  series  was  to  use  a  detonable 
gas  as  the  energy  source  in  the  SALMON- STERLING  cavity  and  to 
fire  the  shots  at  differing  initial  pressures  corresponding  to 
differing  energy  levels.  The  gas  pressure  and  energy  content  of 
the  first  two  shots  fired  in  this  series,  DIODE  TUBE  and  HUMID 
WATER,  were  designed  to  match  the  yield  of  the  STERLING  device 
and  to  provide  a  fully  decoupled  loading  condition.  A  third  shot 
in  the  series,  DINAR  COIN,  was  to  be  fired  with  some  three  times 
the  energy  content  of  the  first  two  shots  and  thus  to  provide  a 
partially  decoupled  loading  condition.  Unfortunately,  technical 
difficulties  leading  to  misfires  and  loss  of  data  on  both  the 
DIODE  TUBE  and  HUMID  WATER  shots  resulted  in  the  cancellation  of 
this  series . 


In  addition  to  the  field  experiments  discussed  above,  at 
least  two  small-scale  experimental  programs  to  study  decoupling 
phenomena  have  been  conducted.  One  of  these  involved  the  detona¬ 
tion  of  small  explosive  charges  in  spherical  cavities  formed  in 
plaster  of  paris  (Reference  5).  Unfortunately,  overburden 
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pressures  were  not  applied  in  conducting  these  experiments; 
consequently,  the  results  are  probably  dominated  by  tensile 
radial  cracking  and  have  the  same  limited  application  as  do  the 
results  of  the  COWBOY  series.  In  addition,  the  diagnostic  tech¬ 
nique  used  in  these  experiments  (electrical  capacitance  probes) 
do  not  provide  good  time-resolved  data  on  a  fundamental  proper¬ 
ty,  such  as  stress  or  particle  velocity,  of  the  generated  shock 
wave . 


A  second  small-scale  experimental  program  has  been  conducted 
at  Physics  International  under  ARPA  support  (Contract  No.  DASA- 
01-70-C-0094)  and  has  resulted  in  the  successful  development  of 
experimental  capabilities  and  techniques  to  study  decoupling 
phenomena  on  a  small  scale.  This  experimental  program  involved 
the  successful  execution  of  several  small-scale  decoupling  experi¬ 
ments  under  conditions  of  high  confining  (overburden)  pressure. 

This  report  discusses  in  detail  the  experimental  apDaratus, 
the  problems  which  had  to  be  overcome,  the  experimental  results, 
and  the  comparison  of  these  results  with  calculated  results.  The 
calculations  described  in  this  report  were  carried  out  to  evaluate 
the  effects  of  yield  model,  equation  of  state  model,  overburden, 
source  amplitude  and  source  frequency  on  decoupling  behavior. 

In  addition  to  the  experimental  results  reported  herein,  a 
series  of  decoupling  calculations  was  performed  using  traditional 
finite  difference  techniques.  A  unique  aspect  of  these  calcula¬ 
tions  was  the  coupling  of  purely  elastic  code,  LEE’K ,  to  the  out¬ 
side  of  the  elastic-plastic  Lagrange  grid.  Because  the  LEEK  code 
is  very  efficient,  the  calculations  could  be  carried  out  to  seis- 
mologically  significant  times  while  retaining  good  resolution  of 
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higher  frequency  phenomena.  In  these  calculations  the  effects  of 
yield  mode 1 ,  equation  of  state  model,  overburden,  source  ampli¬ 
tude  and  frequency  on  decoupling  behavior  were  studied.  Respo 
spectra  analysis  of  the  results  illustrates  how,  in  a  highly 
yielding  and  compactible  media,  energy  in  the  higher  frequency 
modes  is  shifted  into  the  lower  frequency  modes  resulting  in  a 
net  amplification  of  energy  in  the  seismically  important  fre¬ 
quencies  . 
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SECTION  3 

SMALL  SCALE  EXPERIMENTS 


3.1  EXPERIMENTAL  TECHNIQUE 


As  discussed  in  the  preceding  section ,  there  have  been  few 
systematic  experimental  programs  to  study  decoupling  phenomena, 

A  reasonably  broad  experimental  study  of  decoupling  behavior 
incorporating  a  significant  number  of  variables  is  probably 
represented  only  by  the  COWBOY  series  conducted  for  the  ASC  in 
1960  (Reference  3)  .  Although  some  excellent  data  on  decoupling 
behavior  were  obtained  in  SALMON  and  STERLING  nuclear  tests,  the 
high  cost  of  such  nuclear  experiments  precludes  the  systematic 
study  of  decoupling  phenomena. 

If  a  small-scale  experimental  capability  could  be  developed, 
tests  to  evaluate  various  parameters  of  interest  could  be  eco¬ 
nomically  conducted  in  a  laboratory  environment.  Efforts  in  this 
program  were  directed  towards  refining  experimental  techniques 
which  were  partially  developed  for  ARPA  on  Contract  No.  DASA-01- 
6£_C_oi47  during  the  period  August,  1968  through  September,  1969 
(Reference  6).  The  successful  application  of  small-scale  experi¬ 
mental  techniques  to  the  study  of  decoupling  phenomena  would  pro¬ 
vide  data  which  could  be  used  to  critically  evaluate  theoretical 
and  calculational  predictions  at  a  cost  significantly  less  than 
that  required  for  large-scale  field  experiments. 
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In  designing  the  experimental  scale  and  geometry  to  be  used 
in  carrying  out  small-scale  decoupling  experiments,  consideration 
was  given  both  to  the  types  of  large  scale  shots  that  might  be 
meaningfully  modeled  and  to  the  cost  effectiveness  of  various 
experimental  designs.  At  the  time  this  experimental  work  was 
initiated,  the  greatest  amount  of  interest  and  effort  in  decou¬ 
pling  was  directed  towards  the  MIRACLE  PLAY  Series  being  planned 
for  the  SALMON /STERLING  cavity  in  the  Tantum  Salt  Dome.  Conse¬ 
quently,  it  was  decided  to  design  the  small-scale  experiments 
to  m  'del  the  SALMON/STERLING  cavity  in  salt  and  to  use  a  deto- 
nable  gas  (rnethane-oxygen)  as  the  explosive  energy  source.  After 
careful  review,  it  was  concluded  that  a  2-inch-diameter  spherical 
cavity  in  rock  salt  would  be  small  enough  to  be  experimentally 
cost-effective  and  yet  large  enough  to  yield  valid  results. 

In  laying  out  a  small-scale  experimental  study  of  decoupling 
phenomena,  careful  consideration  must  be  given  to  any  scaling 
factors  that  might  influence  the  results.  If  the  peak  and  equi¬ 
librium  pressures  developed  by  an  energy  source  are  not  dependent 
upon  the  size  and  the  total  energy  of  the  source  and  if  the 
response  of  the  material  around  the  source  is  time  independent, 
then  times  and  displacements  will  scale  directly  as  lengths  are 
scaled.  Thus  a  2-inch  diameter  spherical  cavity  would  correspond 
to  a  scale  reduction  of  650  from  the  108-foot-diameter  cavity 
created  in  the  Tatum  Salt  Dome  by  the  SALMON  experiment,  and 
times  and  displacements  at  relative  radii  should  scale  by  a 
factor  of  650.  At  any  given  relative  radius,  the  stresses  and 
particle  velocities  would  be  equivalent  to  those  observed  in  a 
full-scale  experiment. 
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Although  linear  scaling  laws  stipulate  that  experimental 
results  will  not  be  a  function  of  size,  it  is  still  desirable 
to  conduct  small-scale  experiments  on  as  large  a  scale  as  practi¬ 
cal.  If  time  dependent  effects  should  perturb  experimental 
scaling,  the  larger  experiments  would  be  the  more  valid  and 
informative . 

One  of  the  most  desirable  and  most  difficult  parameters  to 
simulate  is  overburden  pressure.  The  only  way  in  which  over¬ 
burden  pressure  could  be  simulated  in  small-scale  decoupling 
experiments  would  be  by  applying  an  appropriate  fluid  pressure 
to  the  exterior  surface  of  the  test  block  containing  the  scaled 
decoupling  cavity.  Either  a  liquid  or  a  gas  may  be  used  to 
provide  this  pressure,  but  in  either  case  an  appropriate  pressure 
vessel  to  contain  the  fluid  must  be  designed  and  the  test  sample 
must  be  properly  encapsulated  so  that  the  fluid  does  not  pene¬ 
trate  the  interstices  of  the  sample  and  the  decoupling  cavity 
within  it.  Since  it  is  desirable  to  protect  the  pressure  vessel 
from  large  shocks  and  overpressures  that  might  be  generated  by 
detonations  within  the  test  decoupling  cavity,  confining  pressures 
are  most  logically  provided  by  a  highly  compressible  inert  gas 
rather  than  a  fluid  such  as  kerosene  or  hydraulic  oil.  Dry 
nitrogen,  which  is  commercially  available  in  prepressurized 
bottles  at  pressures  up  to  6,000  psi,  was  selected  as  the  ideal 
confining  medium. 

A  pressure  vessel  that  had  been  designed  and  built  in  the 
portion  of  this  program  carried  out  on  Contract  No.  DASA  01-68- 
C-0147  was  quite  adequate  for  continuing  the  experimental  studies. 
This  .pressure  vessel,  shown  schematically  in  Figure  1,  is  con¬ 
structed  of  two  pressed  steel  hemispheres  held  together  by  30 
high-scrength  bolts  through  flanges  welded  to  the  hemispheres. 
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Figure  1  Cross  section  of  spherical  pressure  chamber 

containing  test  sample  with  a  2-inch-diameter 
cavity  for  decoupling  detonations. 
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The  most  attractive  feature  of  this  design  is  its  large  internal 
volume  of  spherical  geometry  with  a  diameter  of  22  inches.  This 
internal  geometry  would  contain  a  reasonably  large  equi-dimen- 
sional  test  block,  which  would  provide  the  optimum  in  recording 
times  for  any  instrumentation  used  in  the  small-scale  decoupling 
experiments . 

An  undesirable  feature  of  the  pressure  vessel  is  the  30 
high-strength  bolts  used  to  attach  the  two  hemispheres  together. 
Even  when  optimum  bolt  spacing,  diameter,  and  strength  are  used, 
this  closure  method  greatly  limits  the  maximum  pressures  within 
a  chamber  of  this  design.  The  3000  psi  maximum  pressure  for 
which  this  vessel  was  designed  can  be  obtained  only  by  allowing 
for  a  very  small  (1.1)  safety  factor  in  the  failure  of  the  bolts. 
This  pressure  vessel  could  be  used  at  these  high  pressures  only 
by  employing  remote  operating  procedures  with  the  chamber  c  the 
firing  table  at  Physics  International's  Tracy  Test  Site.  In 
addition,  the  time  to  install  and  torque  up  the  30  bolts  adds 
considerably  to  the  time  required  to  conduct  a  small-scale  de¬ 
coupling  experiment.  A  cylindrical  pressure  vessel  with  threaded 
plug  closures  at  each  end  would  be  a  preferable  design  for  use 
in  any  future  experiments . 


Through  a  series  of  trial  and  error  experiments,  it  was  de¬ 
termined  that  the  high  pressure  nitrogen  used  to  provide  confin¬ 
ing  pressures  could  be  exclude'5  from  the  interstices  of  the  test 
sample  and  from  the  small  scale  decoupling  cavity  by  employing 
concentric  layers  of  two  encapsulating  materials  on  the  surface 
of  the  test  block.  To  provide  strength  and  integrity  to  the 
test  block,  which  was  composed  of  several  carefully  mated  pieces 
of  salt,  a  hard  high-strength  epoxy,  Furane  Epibond  104  with  946 
hardener,  was  used  for  the  first  encapsulating  layer.  To  ensure 
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that  the  high  pressure  nitrogen  could  not  penetrate  through  the 
somewhat  brittle  Epibond  104  layer,  a  second  layer  of  Uralane 
8615  was  applied  after  the  Epibond  layer  had  fully  cured.  Only 
a  few  minor  problems  with  leaks  were  encountered  in  utilizing 
this  two- layer  encapsulating  technique. 

A  photograph  of  the  lower  hemisphere  of  the  high  pressure 
chamber  containing  a  completely  encapsulated  test  block  of  salt 
is  shown  in  Figure  2.  The  tube  coming  out  of  the  top  of  the 
salt  block  is  the  detonable  gas  fill  line.  The  three  white  wires 
on  the  left  side  of  the  salt  block  are  coaxial  strain  gage  leads, 
which  connect  to  the  three  high-pressure  electrical  lead-throughs 
seen  passing  through  the  flange  of  the  lower  hemisphere  on  the 
left  hand  side.  The  steel  tube  seen  exiting  from  the  flange  of 
the  lower  hemisphere  closest  to  the  camera  is  for  the  emplacement 
of  appropriate  detonators  into  the  center  of  the  test  decoupling 
cavity.  Details  of  the  detonable  gas,  strain  gage  diagnostic, 
and  detonator  systems  are  given  below.  After  the  test  block  has 
been  completely  connected  for  the  purposes  of  conducting  a  small- 
scale  decoupling  experiment,  the  upper  hemisphere  is  lowered  into 
place  and  bolted  down  with  the  30  high-strength  bolts  placed 
through  the  holes  shown  on  the  upper  surface  of  the  flange  in 

Figure  2 . 

A  photograph  of  the  assembled  charter  ready  for  the  applica¬ 
tion  of  high  confining  pressure  nitrogen  and  the  execution  of 
the  small-scale  decoupling  experiment  is  shown  in  Figure  3.  The 
smaller  diameter  steel  tubing  seen  running  between  the  pressure 
charter  and  the  remote  valve  system  is  for  the  high  pressure 
nitrogen,  while  the  larger  diameter  tubing  is  for  filling  the 
test  cavity  with  a  detonable  gas  mixture.  The  box  on  the  groun 
under  the  pressure  tubes  contains  the  strain  gage  power  supplies 
and  amplifiers  used  in  measuring  radial  displacements  during  an 

experiment. 
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Figure  3 


Closed  and  bolted  pressure  chamber  ready  for 
pressurization  by  high-pressure  nitrogen  and 
execution  of  a  small-scale  experiment. 
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The  design  of  small-scale  experiments  to  model  the  D/ODE 
TUBE,  HUMID  WATER  and  DINAR  COIN  experiments  of  the  MIRACLE  PLAY 
series  required  that  a  system  for  providing  a  methane-oxygen  mix 
to  the  scaled  decoupling  cavity  be  designed  and  built.  Since 
this  system  had  to  function  only  at  relatively  low  pressures 
(less  than  1000  psi  maximum),  the  design  and  fabrication  of  the 
detonable  gas  fill  system  was  relatively  simple  and  straight¬ 
forward.  To  ensure  that  the  system  could  be  operated  with  a 
high  degree  of  safety,  it  was  necessary  to  use  specially  designed 
explosion-proof  valves-  to  construct  a  remote  fill,  bleed  and 
abort  system  to  operate  on  the  firing  table  at  the  Tracy  Test 
Site,  and  to  provide  valving  to  ensure  that  a  detonation  or 
deflagration  initiated  in  the  test  cavity  could  not  propagate 
back  down  the  lines  of  the  gas-fill  system.  The  remotely  opera¬ 
ted  fill,  bleed,  and  abort  valves  and  an  80 , 000-psi-rated  Auto¬ 
clave  isolation  valve  to  preclude  detonations  from  propagating 
back  into  the  gas-fill  system  are  contained  in  the  steel  box 
seen  partially  on  the  right  of  Figure  3.  With  the  exception  of 
some  leakage  past  the  explosion  proof  valves  at  the  higher  firing 
pressures,  no  difficulties  were  encountered  with  the  detonable 
gas-fill  system. 

At  the  initiation  of  this  research  program  it  was  anticipated 
that  the  transition  from  a  deflagration  to  a  detonation  in  the 
detonable  gas  would  not  naturally  occur  within  the  1-inch  radius 
of  the  scaled  cavity.  It  was  initially  thought,  however,  that 
a  detonation  could  be  readily  initiated  with  a  high  energy  spark 
gap  placed  at  the  center  of  the  cavity.  The  evaluation  of  this 
initiation  system  conducted  on  Contract  DASA  01-68-C-0147 
definitely  indicated  that  detonations  were  not  achieved  (Refer¬ 
ence  6)  and  that  the  energy  of  the  detonable  gas  mix  was  released 
through  a  slowly  propagating  def lagrauion .  Consequently,  efforts 
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to  design  a  high  explosive  detonator  were  bogun  on  Contract 
DASA-01-6  8-C-0147 .  The  design  and  testing  of  a  suitable  deto¬ 
nator  for  these  small-scale  experiments  were  continued  into  the 
current  contract.  A  detonator  that  was  designed  and  success¬ 
fully  tested  is  shown  in  Figure  4.  This  detonator,  containing 
approximately  40  mg  of  PETN  pressed  around  a  centrally  located 
bridgewire  and  encased  in  a  gelatine  capsule,  is  an  improved 
modification  of  the  detonator  originally  developed  on  Contract 
DASA  01-68-C-0147.  This  newer  detonator  geometry  would  assure 
the  generation  of  a  spherically  divergent  detonation  wave,  while 
the  original  detonator  with  an  asymmetrically  located  bridge- 
wire  would  not.  The  2.32  x  109  ergs  energy  contained  in  40  mg 
of  PETN  is  less  than  4  percent  of  the  energy  of  1.06  gm  of  a 
detonable  methane /oxygen  mix  at  220  psi  in  a  2-inch-diameter 
cavity. 

Although  the  PETN  bridgewire  detonators  were  found  to  pro¬ 
vide  very  reliable  detonations  in  me thane -oxygen  gas  mixtures 
ranging  in  pressure  from  110  to  660  psi,  it  was  found  that  they 
could  not  be  used  inside  the  large  spherical  pressure  vessel. 

The  electrical  noise  from  this  detonator  had  been  tolerable  in 
all  tests  conducted  outside  the  pressure  vessel,  but  very  severe 
noise  problems  were  encountered  when  the  detonator  was  fired 
within  the  vessel.  Large  amplitude  electrical  transients  asso¬ 
ciated  with  the  capacitor  discharge  system  used  to  fire  the 
detonator  lasted  for  periods  of  over  100  Msec  when  the  detonator 
was  fired  inside  the  pressure  vessel.  Since  the  shock  wave 
generated  by  the  detonated  gas  will  arrive  at  the  various  gage 
locations  in  the  salt  block,  some  20  Msec  after  the  detonator  is 
fired  the  electrical  noise  would  completely  override  any  gage 
signal.  This  problem  was  readily  resolved  by  placing  the  bridge- 
wire  portion  of  the  detonation  system  outside  the  pressure  vessel 
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Figure  4  Geometry  for  measurement  of  radial  displacements 
with  tangential  strain  gages. 


and  using  one  grain  mild  detonating  fuse  (MDF)  to  carry  the 
detonation  into  the  test  cavity  inside  the  salt  block.  In  order 
to  protect  the  MDF  where  it  passes  through  the  high  confining 
pressure  nitrogen  and  to  reduce  the  mechanical  signal  imparted 
to  the  salt  block  by  the  MDF ,  the  MDF  was  encased  in  adiprene 
inside  a  1/8-inch  o.d.  by  1/16-inch  i.d.  stainless  steel  tube. 

A  40-mg  PSTN  capsule  identical  to  the  one  used  with  the  bridge- 
wire  illustrated  in  Figure  5  was  placed  on  the  end  of  the  MDF 
at  the  center  of  the  test  cavity. 

The  MDF-ini tiated  detonators  were  also  found  to  provide 
very  reliable  initiation  of  detonations  in  methane-oxygen  mix¬ 
tures  with  ratios  of  1:1.6  and  1:2  at  pressures  of  110,  220,  330, 
440,  and  660  psi.  This  pressure  range  completely  covers  the 
initial  gas  pressures  that  would  be  required  for  conducting  both 
underdriven  and  overdriven  small-scale  decoupling  experiments. 

In  utilizing  the  MDF- initiated  detonator  in  carrying  out  the 
small-scale  decoupling  experiments,  the  only  problems  encountered 
were  associated  with  electrical  noise  generated,  as  a  result  of 
ground  loops  created  by  the  inadvertent  grounding  of  the  steel 
tube  containing  the  MDF. 

Several  diagnostic  systems  were  considered  and  carefully 
evaluated  for  possible  use  in  measuring  and  evaluating  decoupling 
effectiveness  of  spherical  cavities  in  salt.  Physics  Interna¬ 
tional  has  had  considerable  experience  in  piezoresistive  and 
piezoelectric  pressure  gages,  electromagnetic  velocity  gages, 
and  strain  gages  in  making  dynamic  measurements  of  material 
response.  Because  they  cannot  be  utilized  as  an  in-material 
gage,  the  piezoelectric  type  gages  were  not  considered.  Piezo¬ 
resistive  gages  such  as  manganin  and  ytterbium  can  be  very 
effectively  used  for  making  pressure  measun  r.ents  at  high  stress 
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Figure  5  PETN-bridgewire  detonator  used  to  initiate  detonations 
in  a  2-inch-diamet€’r  cavity  filled  with  methane  and 
oxygen . 
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levels  (typically  greater  than  1  kbar)  .  Because  the  stress 
levels  anticipated  in  the  small-scale  decoupling  experiments 
would,  in  many  cases,  be  much  less  than  1  kbar,  it  was  not  e 
that  this  type  of  gage  could  be  satisfactorily  used.  At  the 
time  that  a  diagnostics  system  had  to  be  selecteu  or  these 
experiments,  Physics  International's  experience  with  the  electro¬ 
magnetic  velocity  gage  had  not  been  sufficient  to  enable  us  to 
completely  evaluate  its  potential  application  to  this  program. 
Consequently,  dynamic  strain  gages  were  selected  as  the  mos 
promising  diagnostic  system  for  the  experiments. 

In  conducting  the  small-scale  decoupling  experiments,  radial 
displacement  versus  time  histories  could  be  obtained  by  placing 
conventional  dynamic  strain  gages  on  planes  in  the  material  to 
be  studied  at  appropriate  radii  from  the  explosive  cavity, 
illustrated  in  Figure  4,  the  planes  on  which  the  strain  gages 
are  placed  are  normal  to  a  radial  direction  from  the  cavity  and 
the  strain  gages  will  in  effect  measure  the  tangential  or  hoop 
strains  at  that  radius.  If  a  strain  gage  can  make  a  good  measure¬ 
ment  of  the  dynamic  hoop  strains,  then  this  may  be  readily  re¬ 
duced  to  a  radial  displacement  history  given  the  appropriate 
geometrical  considerations.  The  equation  for  reducing  e 
tangential  or  hoop  strains,  cfc.  as  measured  by  a  strain  gage 

l  *  p  nrcurrinq  at  that  radius,  R,  is. 
to  the  radial  displacement,  AR,  occurring  a 

AR  =  R  x  e. 


inclusion  of  the  strain  gage  factor,  G,  and  the  voltage,  V, 
placed  across  the  gage  yields: 


AR  = 


R  x  AV_ 
G  x  V 


where  AV  is  the  change  in  voltage  across  the  gage  caused i  by f  strait n 
induced  changes  in  resistance.  This  equation  was  used  to 
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the  strain  gaqe  voltage  output  obtained  during  an  experiment  to 
radial  displacements.  Note  that,  as  a  strain  gage  is  at  a 
greater  radial  distance  from  a  test  cavity,  the  strain  seen  by 
that  gage  will  go  down  while  the  corresponding  amplification 
factor,  R,  will  increase.  Consequently,  strain  gages  at  large 
relative  radii  to  the  cavity  would  be  expected  to  have  unfavor¬ 
able  signal-to-noise  ratios. 

Strain  gage  amplifiers  with  both  low  and  high  frequency 
responses  were  developed  and  used  to  monitor  strain  gage  output 
on  the  experiments  conducted  in  this  program.  An  oscilloscope 
record  showing  the  response  of  a  single  strain  gage  as  seen  by 
be  th  a  high  frequency  and  a  low  frequency  amplifier  is  shown  in 
Figure  6.  The  high  frequency  amplifier  does  not  seem  to  be 
if fected  by  normal  levels  of  electrical  noise  and  does  show  much 
more  diagnostic  wave  structure  than  does  the  low  frequency  ampli 
fier.  The  high  frequency  amplifier  has  a  maximum  frequency 
response  of  approximately  3  megacycle?  while  the  low  frequency 
amplifier  has  a  maximum  frequency  response  of  approximately 
100  kilocycles.  A  simple  computer  program  was  developed  to 
reduce  the  voltage  output  signal  of  the  strain  gage  amplifiers, 
as  recorded  on  an  oscilloscope  (Figure  6) ,  to  the  displacement 
time  histories  at  that  strain  gage  locations. 

Since  the  inception  of  this  experimental  program,  Physics 
International  has  gained  considerably  more  experience  with 
electromagnetic  velocity  gages  (References  10  and  11)  and  has 
carefully  evaluated  the  possible  applicability  of  this  type  of 
gage  to  making  measurements  of  small-scale  decoupling  experi¬ 
ments.  Although  the  strain  gage  diagnostic  system  did  provide 
good  measurements  of  radial  displacement,  it  would  now  be 
recommended  that  the  electromagnetic  velocity  gages  be  used  m 
any  future  small-scale  decoupling  experiments. 
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Figure  6  Typical  oscilloscope  record  of  strain-gauge  response 
showing  difference  between  low-frequency  and  high- 
frequency  amplifier  output. 
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3.2  RESULTS 


There  are  numerous  parameters  that  might  be  beneficially 
studied  within  the  scope  of  an  experimental  program  such  as  con¬ 
ducted  on  contract  DASA  01-70-C-0094 .  Such  parameters  could 
include  size  and  shape  of  the  decoupling  cavity,  the  explosive 
energy  released  in  that  cavity,  the  strength  and  compacts  ■’  char 
acteristics  of  the  material  forming  the  cavity ,  and  the  over¬ 
burden  pressure  acting  at  the  time  a  decoupling  test  was  con¬ 
ducted.  Since  it  was  not  possible  to  incorporate  all  these 
parameters  within  the  scope  of  the  funded  program,  it  was  neces¬ 
sary  to  decide  which  parameters  could  be  most  meaningfully 
investigated.  Since  the  small-scale  experiments  were  primarily 
designed  to  model  shots  of  the  MIRACLE  PLAY  series,  it  was  felt 
that  the  initial  experiments  should  be  directed  towards  studying 
the  effects  of  increasing  the  explosive  energy  in  a  fixed-size 
spherical  cavity  in  rock  salt.  Once  the  decoupling  effectiveness 
of  the  spherical  cavities  modeling  the  MIRACLE  PLAY  series  had 
been  established,  it  would  then  be  most  desirable  to  study  the 
possible  effects  of  variations  in  material  properties.  Because 
work  hardening  of  the  salt  around  the  SALMON/STERLING  cavity  in 
the  Tatum  salt  dome  could  affect  the  results  of  a  decoupling 
experiment,  it  was  felt  that  a  highly  work-hardened  salt  might  be 
a  very  desirable  medium  in  which  to  conduct  subsequent  shots. 


Nine  small-scale  decoupling  experiments  were  successfully 
conducted  on  the  program  during  the  period  September  28,  1970  to 
December  22,  1970.  All  of  these  experiments  were  carried  out 
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using  a  detonable  me thane -oxygen  gas  mixture  in  a  2-inch-diameter 
spherical  cavity  in  salt.  The  initial  partial  pressures  of 
methane  and  oxygen,  their  ratio,  and  the  total  detonable  gas 
pressure  utilized  in  each  of  the  nine  experiments  are  given  in 
Table  1.  The  first  five  experiments  were  conducted  in  a  coarse¬ 
grained  natural  rock  salt  from  Texas  while  the  last  four  experi¬ 
ments  were  carried  out  in  a  cavity  formed  in  an  artifically 
pressed  salt.  All  of  the  experiments  were  conducted  with  a  con¬ 
fining  pressure  of  2500  psi  designed  to  simulate  the  2720-foot 
depth  of  the  SALMON/ STERLING  cavity. 

Except  for  the  first  and  last  experiments,  excellent  strain- 
gage  records  were  obtained  for  all  strain  gages  (2  or  3)  and  all 
amplifiers  (3  or  4)  used  in  each  experiment.  In  both  the  first 
and  last  experiments,  very  high  level  electrical  noise,  probably 
caused  by  current  ground  loops  formed  when  the  steel  tube,  con¬ 
taining  the  MDF  used  to  fire  the  PETN  capsule  inside  the  cavity, 
was  inadvertently  grounded  to  the  steel  pressure  chamber.  When 
this  grounding  occurred,  current  from  the  capacitor  discharge 
unit  used  to  fire  the  shot  could  travel  down  the  steel  tube , 
through  the  steel  pressure  vessel  and  back  to  ground  through  the 
steel  tubes  of  the  detonable  and  high-pressure-gas  systems. 

Basic  information  obtained  from  the  strain  gage  diagnostic 
system  used  in  the  nine  experiments  conducted  on  this  program  is 
summarized  in  Table  2.  In  this  table  the  numerical  portion  (1, 

2  or  3)  of  the  strain  gage  designation  refers  to  an  actual  strain 
gage  while  the  letter  portion  (L  or  H)  of  the  designation  per¬ 
tains  to  the  type  of  amplifier  used  to  monitor  that  strain 
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TABLE  1 


SUMMARY  OF  EXPERIMENTAL  CONDITIONS  FOR  SMALL-SCALE 
DECOUPLING  EXPERIMENTS 


Shot 

Material 

Firing 
Pressure 
(psi)  (bars) 

240-1 

Natural 

Rock 

Salt 

220 

15.2 

240-2 

II 

440 

30.3 

240-3C 

II 

220 

15.2 

240-4 

II 

4*0 

30.3 

240-5 

•1 

(  60 

45.4 

240-6 

Artifi¬ 

cially 

Pressed 

Salt 

220 

15.2 

240-7 

II 

220 

15.2 

240-8 

II 

440 

30.3 

240-9 

II 

440 

3'  .  3 

Gas 

°2 

(psi) 

Mixture 

ch4 

(psi) 

Ratio 

Confining 

Pressure 

(Psi) _ 

135 

85 

1:1.6 

2500 

270 

170 

1:1.6 

2500 

135 

85 

1:1.6 

2500 

270 

170 

1:1.6 

2500 

405 

255 

1:1.6 

2500 

135 

85 

1:1.6 

2500 

147 

73 

1:2.0 

2500 

293 

147 

1:2.0 

2500 

270 

170 

1:1.6 

2500 

TABLE  2 


SUMMARY  OF  STRAIN  GAGE  RESULTS 
FOR  SMALL-SCALE  DECOUPLING  EXPERIMENTS 


Shot 

Strain  Radius 
Gaqe  (inch) 

Peak 

Displacement 

(cm) 

Arrival  Time 
From  Zero 
Time  (usec) 

Peak 

Arrival  Time 
(Usee) 

From 

Arrival 

Time 

From 

Zero  Tinu 

240 

-1 

1L 

2.960 

0.60  x 

10-4 

104 

3.2 

10  7.2 

2L 

2.913 

0.71  x 

10  4 

104 

3.3 

10  7. 3 

2  40 

-2 

1L 

2.960 

2.61  x 

10-4 

105 

6.2 

111.2 

2L 

2.913 

2.55  x 

10  4 

104 

3.3 

10  7. 3 

240 

-3C 

1L 

3.05 

5.51  x 

10-4 

10  8 

18.  7 

126.7 

2L 

3.12 

5.49  x 

10-4 

111 

12.5 

12  3.5 

2H 

3.12 

7.16  x 

10-4 

10  8 

7.6 

115.6 

1 

3H 

4.60 

3.07  x 

10  4 

114 

5.3 

119.3 

240' 

-4 

1L 

2L 

3.05 

3.12 

11.76  x 
14.07  x 

10-4 

10-4 

116 

118 

18.9 

14.4 

134.9 
132 . 4 

2H 

3.12 

15.93  x 

10-4 

115 

8.8 

12  3.8 

3H 

4.60 

5.31  x 

10  4 

124 

5.5 

129.5 

240- 

-5 

1L 

3.05 

16.2  3  x 

10-4 

120 

12.0 

126.0 

2L 

3.12 

16.89  x 

10-4 

120 

14.8 

132.9  1 

2H 

3.12 

17.98  x 

10  _4 

119 

12.6 

131.0 

3H 

4.60 

5.87  x 

10  4 

126 

7.0 

133.0 

240- 

-6 

1H 

3.008 

5.33  x 

10-4 

100 

8.0 

10  8.0 

2H 

4.041 

6.24  x 

10-4 

10  7 

8.5 

115.5 

2L 

4.041 

3.81  x 

10-4 

108 

14.5 

122.5 

3H 

6.0  75 

4 . 15  x 

10  4 

120 

8.5 

128.5 

240- 

•7 

1H 

3.008 

7.36  x 

10-4 

105 

8.1 

113.1  | 

2H 

4.041 

7.58  x 

10-4 

112 

7.7 

119.7 

2L 

4.041 

4.12  x 

10-4 

112 

12.0 

12  4.0 

3H 

6.075 

5.75  x 

10  4 

129 

8.5 

137.5 

I 

240- 

8 

1H 

3.008 

14.10  x 

iU-4 

102 

9.2 

111.2 

2H 

4.041 

12 . 46  x 

10-4 

108 

8.3 

116.3 

2L 

4.041 

8.00  x 

i0_4 

109 

14.5 

12  3.5 

3H 

6.075 

9.20  x 

10  4 

120 

8.6 

12  8.6 

2  40- 

9 

1H 

3.008 

12.30  x 

10-4 

120 

8.0 

12  8.0 

2H 

4.041 

12.83  x 

10_4 

128 

12.0 

140.0 

3H 

6.075 

9.50  x 

10  4 

134 

7.0 

141.0 

30 
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gage.  The  letter  "L"  indicates  that  a  low-frequency,  100-kHz-- 
peak- response  strain  gage  amplifier  was  used  whereas  the  letter 
"H"  indicates  a  high-frequency,  3-MHz  strain  gage  amplifier. 

Note  that  in  nearly  all  shots  one  .strain  gage  was  monitored  with 
both  a  low-  and  a  high-frequency  amplifier  so  that  relative 
response  data  for  these  two  amplifying  systems  Could  be  obtained. 
The  time  of  arrival  of  the  initial  shock  wave,'  the  peak  dis¬ 
placement  observed  by  each  strain  gage,  and  the  arrival  time  of 
that  peak  displacement  are  indicated  in  Table  2  also.  A  dif¬ 
ference  in  amplitude  and  arrival  times  of  the  peaks  for  the  low-  , 
and  high-frequency  amplifier  records  of  a  single  strain  gage  is 
due  to  the  differences,  in  responses  of  these  amplifiers  (see, 
for  example.  Figure  6) .  The  arrival  time  for  the  initial  shock 
at  each  strain  gage  includes  time  for  the  detonation  of,the,MDF, 
detonation  of  the  detonable  gas  across  the  cavity,  and  the 
propagation  of  the  shock  wave  from  the  cavity  wall  to  the 
referenced  strain  gage.  The  differences  in  arrival  times  for 
strain  gages  at  different  radii  did  provide  data  which  icould  be 
used  to  calculate  the  shock  velocity  in  the  salt  material  being 

studied.  1  ; 

•  .  •  1 
;  ,  I 

The  first  small-scale  decoupling  experiment  (240-1)  was 
conducted  on  Septepnber  28,  19  70.  A  220-psi  mixture  of  methane 
and  oxygen  in  a  ratio  of  1  to  1.6  was  used  as  the  detonable  gas 
energy  source  in  this  experiment.  With  the  exception  of  some 
adverse  electrical  noise  associated  with  ground  loops  formed  via 
a  shorted  MDF  tube,  every  aspect  of  the  system  appeared  to  work 
quite  reliably  during  the  execution  of  this  first  test  shot. 
Despite  this  noise,  interp'retable  records  were  obtained  from  both 
of  the  strain  gages  used  in  this  shot.  In  this  experiment  only 
the  low-frequency  amplifiers  were  utilized.  Tl^e  reduced  dis¬ 
placement  time  histories  recorded  by  the  two  strain  gages  are 
shown  in  Figure  7.  A  sharp  increase  in  displacement  shown  at  tfie 
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end  of  each  record  is  due  to  the  arrival  of  reflected  waves  at 
the  gage  and  marks  the  end  of  valid  gage  recording.  The 
significance  of  these  records  will  be  discussed  in  raore  detail 
below  along  with  a  discussion  of  the  records  obtained  in  Shot 
240-2,  which  are  also  shown  in  Figure  7, 

One  of  the  most  rewarding  results  of  Shot  240-1  wa^  the 
post-shot  condition  of  the  shot  block  cavity  and  jacketing 
material.  After  the  shot  was  fired,  the  detonable-gas  fill  line 
was  vented  and  the  pressure  in  the  cavity  monitored.  There  was 
no  evidence  of  any  leakage  of  the  high  pressure  confining  gas 
through  the  jacketing  materials  into  the  salt  cavity,  indicating 
that  damage,  if  any,  to  the  salt  block  was  minor.  After  the 
salt  block  was  removed  from  the  high-pressure  chamber,  it  was 
disassembled  along  the  central  parting  plane  so  that  the  cavity 
and  the  gas  fill  and  detonator  lines  could  be  inspected.  A  photo¬ 
graph  of  the  post  shot  cavity  after  Shot  240-1  is  shown  in 
Figure  8.  A  lack  of  damage  to  the  salt  immediately  adjacent  to 
the  cavity  and  the  small  extent  to  which  the  detonation  products 
escaped  along  the  parting  plane  both  attest  to  the  fact  that  the 
cavity  was  underdriven  by  the  detonation. 

The  excellent  condition  of  the  post-shot  salt  block  and  the 
long  time  response  of  the  strain  gages  (including  the  recording 
of  reflected  waves)  indicated  that  the  salt  block  could  be  re¬ 
assembled  for  a  second  shot.  This  reuse  of  the  salt  block  for 
sequential  shots  is  especially  valuable  in  terms  of  gaining  high 
quality  data.  As  discussed  above,  large  grain  size  of  the  salt 
and  the  strain-gage  coupling  to  the  salt  can  have  an  effect  upon 
actual  strain-gage  response.  The  difference  in  the  response  be¬ 
tween  two  strain  gages  with  essentially  identical  locations  in  the 
salt  block  used  for  Shot  240-1  indicated  that  these  effects  could 
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be  fairly  serious.  If  a  single  strain  gage  installation  could 
be  used  for  several  shots,  however,  then  the  effects  of  large 
grain  size  and  strain-gage  coupling  would  be  eliminated  as  far 
as  the  relative  response  of  that  strain  gage  was  concerned. 
Because  it  was  not  necessary  to  change  or  modify  the  strain  gages 
used  in  Shot  240-1  in  refabricating  the  salt  block  for  Shot 
240-2,  a  very  valid  comparison  of  the  strain-gage  records  ob¬ 
tained  in  the  two  shots  may  be  made. 

The  displacement- time  histories  measured  by  the  two  strain- 
gages  in  Shot  240-2  are  shown  in  Figure  7  with  the  records  for 
Shot  240-1. 

Shot  240-2  was  fired  using  a  440-psi  mixture  of  CH4  and  C>2 
in  the  ratio  of  1  to  1.6.  Since  this  pressure  is  exactly  twice 
that  used  for  Shot  240-1,  this  shot  had  essentially  twice  the 
explosive  energy  of  the  earlier  shot.  Displacement  profiles 
from  these  two  shots,  as  illustrated  in  Figure  7,  indicate  that 
displacements  increased  by  more  than  a  factor  of  three  between 
the  two  shots.  This  effect,  or  course,  is  exactly  what  one 
would  expect  in  going  from  a  fully  decoupled  shot  to  a  partially 
decoupled  shot,  as  was  the  case  for  Shot  240-2. 

The  excellent  condition  of  the  salt  block  after  Shots  240-1 
and  240-2  lead  us  to  positively  conclude  that  a  single  salt 
block  could  be  used  quite  successfully  for  several  shots.  Con¬ 
sequently,  we  redesigned  our  experimental  arrangement  so  that, 
once  a  salt  block  was  fitted  with  strain  gages  and  sealed  to 
withstand  the  high  confining  pressures,  it  could  be  refitted 
with  the  necessary  detonator  without  having  to  be  removed  from 
the  large  hemispherical  pressure  chamber.  The  modified  con¬ 
figuration  involved  using  a  large -diameter  (3/8-inch  o.d.  x 
1/4-inch  i.d.)  steel  tube  to  go  from  the  outside  pressure  chamber 
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to  inside  the  salt  cavity.  The  salt  blocks  used  for  the  second 
and  third  series  of  shots  employed  this  new  configuration.  The 
salt  block  used  for  the  second  series  of  experiments  was  also 
manufactured  of  the  same  coarse-grained  Texas  salt  dome  material 
as  was  used  for  the  first  two  shots.  As  in  the  configuration 
for  the  first  salt  block  shots,  two  strain  gages  were  located  at 
approximately  three  cavity  radii  from  the  center  of  the  shot.  A 
third  strain  gage  was  located  at  approximately  4.6  cavity  radii 
from  the  center  of  the  cavity.  Other  aspects  of  the  configura¬ 
tion  of  this  salt  block  will  be  discussed  in  greater  detail 
later. 

In  firing  Shots  240-3,  240-4,  and  240-5,  efforts  were  made 
to  obtain  three  distinctly  different  preshot  cavity  pressures. 

It  was  decided  that  preshot  pressures  of  220  psi ,  440  psi,  and 
660  psi  would  yield  the  most  valuable  data.  Excellent  strain- 
gage  records  were  obtained  from  all  three  strain  gages  for  all 
three  of  these  shots.  Difficulties  with  electrical  noise 
required  that  Shot  240-3  (at  220  psi)  be  refired  several  times, 
until  good  records  w eyp  obtained. 

The  first  successful  data  shot  at  this  pressure  has  been 
designated  as  240-3C.  A  second  shot  in  this  series,  240-4,  was 
successfully  fired  without  any  difficulties  at  440  psi.  Some 
difficulty  was  encountered  with  the  detonatable-gas,  mixed 
pressurization  system  in  attempting  to  fire  the  third  shot  in 
the  series  (Shot  240-5).  A  serious  leak  developed  in  the 
pressurization  system  at  this  higher  pressure.  We  had  difficulty 
maintaining  the  660-psi  pressure  and  could  not  be  sure  that  the 
gas  mix  was  exactly  1.6  02  to  1.0  CH4 .  All  of  the  shots  in  this 
serier  were  fired  with  methane-oxygen  in  the  ratio  of  1  to  1.6, 
and  axl  the  shots  had  an  overburden  pressure  of  2500  psi. 
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Strain-gage  amplifiers  with  both  low-frequency  and  high-frequency 
response  were  used  in  this  sho:.  Strain  gage  1  was  at  a  radius 
relative  to  the  cavity  of  3.05  and  was  monitored  with  a  low- 
frequency  amplifier;  strain  gage  2  at  a  relative  radius  of  3.12 
was  monitored  with  both  low-  and  high  frequency  amplifiers;  and 
strain  gage  3  at  a  relative  radius  of  4.60  cavity  radii  was 
monitored  with  only  a  high-frequency-response  amplifier.  The 
relative  response  of  strain  gage  2,  as  seen  by  both  a  low- 
and  high-frequency  amplifier  in  Shot  240-3C,  is  illustrated  in 
Figure  6. 

The  displacement-time  profiles,  as  monitored  by 
strain  gages  1  and  2  for  the  three  successful  data  shots  in 
this  series,  are  shown  in  Figure  9.  All  six  of  these  profiles 
were  reduced  from  low-frequency  amplifier  records, 
displacement-time  profiles  obtained  from  strain  gages  2  and  3 
using  high-frequency  strain-gage  amplifiers  are  illustrated  in 
Figure  10.  Note  that  the  high-frequency  amplifiers  generally 
show  the  same  response  as  the  low-frequency  amplifiers  but  with 
greater  wave  detail.  As  was  the  case  for  the  first  shots  fired 
in  this  program  (240-1  and  240--2)  ,  tha  peak  displacement  measured 
by  the  strain  gages  increases  by  a  factor  greater  than  two  when 
the  explosive  energy  in  the  cavity  is  doubled  (by  increasing 
initial  gas  pressures  from  220  psi  to  440  psi).  The  very  slight 
increase  in  maximum  displacement  obtained  in  Shot  240-5  is 
probably  a  result  of  having  a  poor  detonatable  gas  mix  and  less 
than  the  anticipated  660-psi  preshot  cavity  pressure. 

Comparison  of  the  strain-gage  records  obtained  in  the  first 
series  of  shots  (240-1  and  240-2,  Figure  7)  with  the  strain-gage 
records  obtained  in  the  second  series  of  shots  (240-3,  240-4, 
and  240-5,  Figure  9)  at  an  approximate  radius  of  3  inches  from 
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Figure  10  Reduced  ("isplacement  versus  time  profiles  for  strain 
gages  2  and  3  as  monitored  by  high  frequency 
amplifiers  in  Shots  240-3c,  240-4  ,  and  240-r,. 


the  center  of  the  cavity  indicates  some  possibly  large  dis¬ 
crepancies  between  the  data  obtained  in  the  two  series.  The 
displacement- time  histories  obtained  in  the  first  series  of  shots 
are  both  lower  in  amplitude  and  much  narrower  in  time  than  in 
the  second  series  of  shots.  These  discrepancies  can  be  most 
readily  attributed  to  the  way  in  which  the  salt  blocks  were 
fabricated  for  each  of  the  series.  The  geometry  of  the  salt 
block  configuration  used  in  the  first  series  is  illustrated  in 
Figure  11.  In  this  block  the  parting  plane  used  to  provide 
access  to  the  cavity  and  the  necessary  gas-filling  and  detonating 
equipment  was  uncemented  when  the  salt  block  was  finally  as¬ 
sembled.  In  addition,  the  planes  on  which  the  strain  gage  was 
located  were  oriented  at  45  degrees  to  this  uncemented  parting 
plane.  If  the  pressure  of  the  detonated  gas  in  the  cavity  were 
high  enough  to  create  tensile  hoop  stresses  around  the  cavity, 
the  uncemented  parting  plane  would  separate  with  a  simultaneous 
injection  of  high-pressure  detonated-gas  products  into  the 
parting  plane.  As  may  be  seen  in  Figure  11,  the  signal  generated 
by  the  injection  of  hot  gases  along  the  parting  plane  would 
arrive  at  the  strain-gage  location  very  soon  after  the  signal 
traveling  along  the  shortest  radial  direction.  Since  the 
anomalous  signal  originating  at  the  uncemented  parting  plane 
would  destroy  the  spherical  geometry  of  an  outgoing  shock  wave, 
it  could  cause  both  the  low  amplitude  and  narrow  time  width  of 
the  strain  gage  profiles  obtained  in  the  first  series  of  shots. 

In  contrast,  the  salt  block  assembled  for  the  second  series 
of  shots  had  a  geometry  less  amenable  to  such  anomolous  effects. 
The  overall  geometry  of  the  salt  block  for  this  series  is  shown 
in  Figure  12.  Because  the  large -diameter  detonator  pipe  allowed 
for  the  replacement  of  the  detonator  at  the  center  of  the  cavity 
without  disassembling  the  salt  block,  we  were  able  to  permanently 
cement  the  central  parting  plane  with  a  high-strength  epoxy. 
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Figure  11  Schematic  of  salt  block  used  for  Shots  240-1  and 

240-2,  showing  the  location  of  strain  gages  around 
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Post-shot  evaluation  of  the  salt  block  after  firing  shots 
240-3,  240-4,  and  240-5  indicated  that  this  cemented  assembly 
plane  had  maintained  its  integrity  during  the  shots.  There  was 
very  little  evidence  of  injection  of  hot  gaseous  detonation 
products  along  this  plane  as  compared  to  the  salt  block  used  in 
the  first  series.  In  addition,  the  location  of  strain  gages  in 
this  salt  block  is  preferable  to  those  used  in  the  salt  block 
of  the  first  series.  As  shown  in  Figure  12,  the  strain  gages 
are  all  located  on  planes  parallel  to  the  central  parting  plane. 
If  there  was  any  failure  of  the  cemented  parting  plane  and 
associated  injection  of  high-pressure  gas  into  this  p..ane,  then 
the  anomolous  signal  so  generated  would  have  less  tendency  to 
perturb  the  spherically  divergent  shock  wave  seen  by  the  strain 
gages.  Consequently,  it  is  believed  that  the  larger-amplitude 
and  broader- wave  profiles  observed  in  the  second  series  of  shots 
are  more  representative  of  true  behavior. 

Because  of  some  questions  about  the  effects  of  the  large 
grain  size  of  the  Texas  salt  dome  material  upon  strain  gage 
coupling  and  response  and  in  order  to  investigate  the  effects  of 
work  hardening  upon  decoupling ,  it  was  decided  that  artificially 
pressed  fine-grain  salt  should  be  used  for  the  third  series  of 
small-scale  decoupling  experiments.  Two  large  blocks  of  artifi¬ 
cially  pressed  fine-grain  salt  were  prepared  by  Autoclave 
Engineers  of  Erie,  Pennsylvania.  Each  of  these  blocks  was 
cylindrical  in  shape  with  an  outside  diameter  of  approximately 
9  inches  and  a  length  of  approximately  18  inches.  One  of  these 
blocks  had  a  conical  opening  with  a  2-inch-diameter  hemisphere 
at  its  bottom  and  pressed  into  the  salt  block.  To  eliminate 
machining  of  a  test  cavity,  this  block  was  used  for  the  third 
series  of  experiments. 
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The  overall  geometry  and  the  straii  gage  locations  for  the 
artificially  pressed  salt  block  used  :n  -he  third  series  of  small 
scale  decoupling  experiments  is  shown  schematically  in  Figure  13. 
In  this  block  all  three  strain  gages  were  located  along  a  single 
radial  direction  at  radial  distances  from  the  center  of  the 
cavity  of  approximately  3,  4,  and  6  inches.  A  spherical  decoupl¬ 
ing  cavity  was  created  in  this  block  by  machining  a  tapered  plug 
of  aluminum  to  fit  the  tapered  hole  pressed  into  the  block.  As 
indicated  in  Figure  13,  a  hemispherical  cavity  was  machined  in 
the  end  of  the  tapered  aluminum  plug  and  the  1/4-inch-i .d.  tube 
used  to  provide  detonator  access  to  the  center  of  the  cavity  was 
placed  in  the  center  of  this  plug.  The  detonable-gas  fill  line 
was  inserted  into  the  test  block  through  a  hole  drilled  into  the 
cavity  from  a  cylindrical  outside  surface  of  the  test  block. 


The  excellent  data  obtained  with  the  high-frequency  strain 
gage  amplifiers  in  earlier  shots  indicated  that  this  type  of 
amplifier  should  be  used  on  all  strain  gages  in  this  third  series 
of  small-scale  decoupling  experiments.  Consequently,  the  strain 
gage  amplifier  box  was  modified  to  contain  three  high-frequency 
amplifiers  and  just  one  low-frequency  amplifier.  Strain  gage  2 
was  monitored  with  both  a  high-frequency  and  a  low-frequency 
amplifier  in  experiments  240-5  through  240-9. 


In  order  to  evaluate  the  effect  of  various  gas  mixtures 
upon  the  characteristics  of  detonation  and  the  possible  quality 
of  data  in  both  small  and  large  scale  experiments,  tests 
utilizing  different  gas  mixtures  were  included  in  the  third 
series  of  experiments.  Test  Shots  240-6  and  240-9  used  methane 
and  oxygen  in  the  ratio  of  1:1.6,  as  in  earlier  experiments.  In 
Shots  240-7  and  240-8,  however,  a  methane-oxygen  ratio  of  1:2 
was  used.  All  of  these  experiments  were  conducted  with  a  con¬ 
fining  pressure  of  2500  psi,  as  were  all  preceding  experiments. 
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Displacement-time  histories  obtained  from  the  three 
strain  gages  in  Shots  240-6  and  240-9  are  illustrated  in  Figure 
14.  The  strain  gage  records  obtained  in  Shot  240-6  were 
generally  of  high  quality  and  had  a  character  in  general  agree¬ 
ment  with  the  displacement- time  histories  observed  in  earlier 
experiments.  The  peak  radial  displacement  monitored  by  the 
strain  gages  in  Shot  240-6  fired  with  an  initial  detonable  gas 
pressure  of  220  psi,  was  in  very  close  agreement  with  that 
obtained  in  a  comparable  shot  fired  a  coarse-grained  natural 
rock  salt  (Shot  240-3c) .  Because  of  its  large  radial  distance 
from  the  cavity  (approximately  6.0  inches)  strain  gage  3 
exhibited  a  low  signal -to -noise  ratio  as  discussed  earlier,  and 
the  record  for  this  strain  gage  should  not  be  considered  as 
valid  as  those  for  strain  gage  1  and  2.  In  Shot  240-9,  which 
was  fired  with  an  initial  gas  pressure  of  440  psi,  significant 
amounts  of  electrical  noise  were  observed  on  the  strain  gage 
oscilloscope  records.  This  noise  was  again  believed  to  be  due 
to  the  generation  of  ground  loops  caused  by  shorting  of  the 
steel  MDF  tube  to  the  steel  high-pressure  chamber.  This  elec- 
trical  noise  was  such  that  the  strain  gage  records  for  gages  1 
and  3  are  of  minimal  value.  As  shown  in  Figure  14,  however,  a  good 
record  was  obtained  for  strain  gage  2.  The  peak  radial  displace- 
ment  observed  by  strain  gage  2  in  Shot  240—9  does  not  show  the 
marked  increase  in  decoupling  which  had  been  observed  for  shots 
fired  at  this  pressure  in  the  coarse-grained  salt.  The  signifi¬ 
cance  of  this  result  will  be  discussed  in  more  detail  later. 

Two  experiments,  240-7  and  240-8,  were  conducted  using 
methane  and  oxygon  in  the  ratio  of  1:2.  Excellent  strain  gage 
records  were  obtained  from  all  strain  gages  in  both  of  these 
experiments.  Displacement- time  histories  obtained  from 
the  three  strain  gages  for  the  two  shots  conducted  with  initial 
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detonable  gas  pressures  of  220  and  440  psi  are  illustrated  in 
Figure  15.  The  displacement-time  records  obtained  from  the  shot 
fired  at  220  psi  (shot  240-7)  show  peak  displacements  appreciab’v 
greater  than  had  been  observed  in  earlier  experiments  conducted 
at  this  initial  firing  pressure.  The  greater  peak  displacement 
amplitude  observed  for  a  shot  fired  with  a  1:2  detonable  gas 
ratio  may  indicate  that  this  mix  is  more  nearly  stoichiometric 
than  the  1:1.6  ratio  for  detonation  under  these  initial  condi¬ 
tions  . 

The  most  interesting  result  obtained  from  the  two  experi¬ 
ments,  240-7  and  240-8,  is  the  apparent  lack  of  increased  coupling 
associated  with  an  increase  in  the  initial  detonable  gas  pressure 
to  a  level  which  had  caused  significant  increases  in  coupling  in 
earlier  experiments.  As  indicated  in  Figure  15,  the  peak  radial 
displacements  observed  in  Shot  240-8  are  not  even  double  those 
observed  in  Shot  240-7,  although  the  initial  gas  pressure  is 
doubled.  It  is  hypothesized  that  this  lack  of  increased  coupling 
is  related  to  the  work-hardened  nature  of  the  salt  media  in  which 
the  tests  were  fired. 

Increased  coupling  can  occur  only  if  the  cavity  begins  to 
respond  non-elastically  as  the  energy  released  by  detonations 
within  it  is  increased.  Generally  such  non-elastic  deformations 
can  be  effected  either  by  tensile  radial  cracking  associated  with 
large  hoop  strains  as  occurred  in  the  COWBOY  experiments  or  by 
extensive  plastic  deformation  of  the  material  forming  the  cavity 
as  occurred  in  the  SALMON  shot  and  as  probably  occurred  in  the 
earlier  high-pressure  shots  utilizing  a  coarse-grained  salt  in 
this  experimental  program.  The  high  confining  pressure  (2500 
psi)  used  in  these  experiments  could  serve  to  preclude  tensile 
radial  cracking,  especially  at  the  stress  levels  generated  by  the 
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Figure  15  Reduced  displacement  versus  time  profiles  obtained 
by  strain  gages  at  three  radii  in  Shots  240-7  and 
240-8. 
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detonable  mixtures  used.  Because  the  artificially  pressed  salt 
was  subjected  to  extensive  deformation  and  work-hardening  during 
the  pressing  operation,  it  would  be  expected  that  the  stress 
levels  to  initiate  plastic  flow  in  this  material  would  be 
significantly  higher  than  for  the  naturally  annealed  coarse- 

grained  salt. 
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3.3  DATA  ANALYSIS 


Analysis  and  evaluation  of  the  experimental  data  were  con¬ 
ducted  using  three  different  approaches.  As  indicated  in  the 
preceding  section,  the  peak  radial  displacements  as  measured  by 
strain  gages  at  various  locations  in  the  test  salt  blocks  were 
compared  to  an  correlated  with  the  explosive  energy  in  the  cavity 
and  the  material  in  which  the  test  was  conducted.  In  an  effort 
to  more  quantitatively  define  the  displacement  response  at  the 
various  strain  gage  locations,  all  high-quality  displacement-time 
profiles  were  subjected  to  a  spectral  analysis.  Finally,  the 
displacement-time  data  and  the  associated  response  curves 
obtained  from  the  small-scale  decoupling  experiments  were  appro¬ 
priately  scaled  and  compared  with  the  results  of  one-dimensional 
finite-difference  calculations  made  for  a  broad-base  parameter 
study  of  large-scale  decoupling  phenomena. 

Variations  in  peak  radial  displacement,  as  observed  by 
strain  gage  measurements  at  several  different  radii  in  the  nin 
decoupling  experiments  carried  out,  are  generally  in  good  agree¬ 
ment  with  anticipated  results.  Any  significant  deviations  from 
expected  behavior  can  be  attributed  to  experimental  difficulties, 
such  as  the  poor  detonable  gas  mix  that  probably  existed  in  Shot 
240-5.  Increases  in  peak  radial  displacement  by  a  factor  of  3.5 
in  the  doubling  of  explosive  energy  of  Shot  240-2  as  compared  to 
240-1,  and  a  factor  of  2.5  in  Shot  240-4  as  compared  to  240-3c, 
experimental  verification  of  the  large  scale  cavity 
decoupling  theory.  In  both  of  these  higher-energy  exper!  ents 
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(240-2  and  240-4)  ,  the  naturally  annealed  coarse  grained 
salt  would  be  expected  to  yield  plastically  at  the  higher  stress 
levels  generated  at  the  cavity  wall.  According  to  the  theory  of 
large-cavity  decoupling,  such  yielding  at  the  cavity  wall  would 
increase  the  radius  at  which  all  deformation  became  elastic  and 
would  accordingly  increase  radial  displacement  seen  at  any 
radius  over  that  which  would  be  seen  if  the  wall  of  the  cavity 
were  to  respond  elastically.  It  is  unfortunate  that  the  pressure 
and  mix  of  the  detonable  gas  used  in  Shot  240-5  are  not  better 
known  and  that  more  experiments  in  the  coarse-grained  salt  could 
not  be  conducted  at  higher  firing  pressures. 

Peak  radial  displacement  measurements  obtained  in  the  third 
series  of  shots  (240-G  through  240-9)  indicate  that  a  stronger 
material  can  provide  a  better  decoupling  behavior  independent 
of  the  effects  of  overburden  pressure.  Comparison  of  peak  radial 
displacements  obtained  in  Shot  240-8  with  those  from  240-7 
indicates  that  no  increased  coupling  occurred  with  a  doubling 
of  the  explosive  energy  content  within  the  cavity.  This  effect 
may  be  readily  attributed  to  the  higher  strength  resulting  from 
work-hardening  of  the  salt  associated  with  the  artificial  press¬ 
ing  of  the  test  blocks.  The  one;  good  strain  gage  record  obtained 
in  Shot  240-9  does  indicate  some  increase  in  coupling  associated 
with  a  doubling  of  the  explosive  energy  as  compared  to  Shot  240-6. 
This  increased  coupling  by  the  small  factor  of  2.2  might  be 
attributable  to  a  weakening  of  the  salt  either  by  damage  asso¬ 
ciated  with  the  firing  of  the  three  previous  shots  (240-6  through 
240-8)  or  the  thermal  annealing  of  the  salt  by  the  hot  detonation 
products  of  these  earlier  experiments. 
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In  carrying  out  this  third  series  of  experiments,  it  was 
hoped  that  at  least  one  experiment  could  be  conducted  with  an 
even  higher  energy  source,  such  as  would  be  obtained  with  a 
660-psi  deton able-gas  mix  in  the  test  cavity.  Unfortunately,  a 
severe  fracture  in  the  salt  block  was  formed  during  the  execu¬ 
tion  of  Shot  240-9.  Besides  damaging  the  salt  block  for  possi¬ 
ble  use  in  further  experiments,  this  fracture  also  served  to 
rupture  the  impermeable  coa'ings  around  the  salt  block,  thus 
allowing  the  high-pressure  nitrogen  confining  media  to  penetrate 
the  experimental  cavity.  An  experiment  with  a  higher  energy 
source  would  serve  to  further  confirm  and  possibly  extend  the 
effects  of  high-strength  material  properties  upon  improved 
decoupling  performance. 

Observation  of  the  peak  radial  displacements  observed  by 
various  strain  gages  is  not  necessarily  the  most  valid  way  of 
analyzing  decoupling  phenomena.  The  peak  displacements  observed 
at  any  strain  gage  location  will  be  the  result  of  the  additive 
effect  of  displacements  occurring  over  a  broad  range  of  fre¬ 
quencies.  Since  the  very  high  frequency  portions  of  the  dis¬ 
placement  waves  will  not  be  propagated  for  large  seismic 
distances,  they  are  not  particularly  pertinent  to  the  decoupling 
problem.  Only  lower- frequency  components  of  the  displacement 
history  are  of  particular  interest. 

In  order  to  critically  evaluate  the  frequency  content  of 
the  displacement-time  profiles  obtained  at  various  gage  loca¬ 
tions,  a  spectral-analysis  computer  program  was  modified  and 
used  to  obtain  response  spectra.  Peak  displacement  and  accelera¬ 
tion  response  spectra  and  Fourier  amplitude  spectra  were  obtained 
for  all  strain  gage  records  possessing  enough  quality  to  ot 
amenable  to  the  analysis.  The  response  spectra  represent  the 
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actual  peak  values  of  displacement,  and  acceleration  obtained 
during  the  analysis  of  a  displacement-time  profile.  Because  of 
the  special  interest  in  the  peak  radial  displacement  as  a  func¬ 
tion  of  frequency  in  the  decoupling  problem,  and  because  the  dis¬ 
placement  response  spectra  were  most  indicative  of  differences 
between  various  displacement  in  time  profiles,  nearly  all  the 
analysis  was  directed  towards  the  study  and  evaluation  of  dis¬ 
placement  response  spectra.  A  discussion  of  the  spectral 
analysis  method  is  given  in  Section  4.1.3  and  in  Figure  45. 

In  order  that  the  response  spectra  obtained  from  the  small- 
scale  data  could  be  more  meaningfully  compared  with  the  results 
of  large-scale  parametric  calculations  and  actual  field  measure¬ 
ments  ,  the  experimental  displacement-time  records  were  appro¬ 
priately  scaled.  As  has  been  noted,  the  small-scale  experiments 
were  reduced  by  a  factor  of  650  from  the  MIRACLE  PLAY  field 
experiments.  Both  displacement  and  time  would  have  to  be  multi¬ 
plied  by  650  to  scale  the  small-scale  results  for  comparison  with 
full-scale  calculations  or  field  data.  All  displacement-versus- 
time  data  from  the  small-scale  experiments  were  scaled  by  a 
factor  of  650  before  spectral  analyses  were  performed.  Conse¬ 
quently,  the  peak  displacement  amplitudes  and  the  frequencies 
indicated  in  Figures  16  through  20  are  comparable  to  the  large- 
scale  calculational  results  of  Figure  21. 

Due  to  the  large-amplitude  electrical  noise  observed  on  the 
strain-gage  records  for  Shot  240-1,  these  records  were  not  suit¬ 
able  for  a  meaningful  spectral  analysis.  Consequently,  no  efforts 
were  made  to  conduct  spectral  analyses  of  the  reduced  displace- 
ment-versus-time  profiles  obtained  in  this  experiment.  Spectral 
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analyses  were  conducted  on  the  scaled  displacement  profiles  from 
Shot  240-2.  These  response  spectra  were  in  general  agreement 
with  those  obtained  on  profiles  from  other  shots,  but  since  there 
are  no  spectra  with  which  they  may  be  directly  compared,  they  are 
not  shown  here. 

The  displacement  response  spectra  obtained  by  applying  the 
analysis  to  the  displacement  profiles  obtained  from  strain  gage 
2L  in  Shots  240-3c,  240-4,  and  240-5  are  illustrated  in 
Figure  16.  As  noted  in  Section  4,  the  strain  gage  displacement 
profiles  were  believed  to  be  valid  for  some  25  microseconds  of 
recording  time.  When  scaled  by  a  factor  cf  650  to  correspond  to 
the  dimensions  of  the  SALMON/ STERLING  cavity,  these  recording 
times  would  correspond  to  a  record  of  some  17  milliseconds  dura¬ 
tion.  Since  wavelengths  with  periods  more  than  twice  as  long  as 
the  maximum  recording  time  would  not  be  adequately  resolved,  the 
values  of  peak  spectral  displacement  for  frequencies  below  30  Hz 
could  be  higher  than  calculated.  The  most  significant  feature  of 
the  response  spectra  illustrated  in  Figure  16  is  the  factor  of 
2.57  difference  at  10  to  70  Hz  between  the  spectra  obtained  from 
strain  gage  2L  in  Shots  240-3c  and  240-4,  Since  the  detonable 
gas  energy  in  the  cavity  was  only  twice  as  high  in  240-4  as  it 
was  in  240-3c,  these  response  spectra  confirm  the  increased 
coupling  associated  with  the  higher-energy  source  used  in  Shot 
240-4.  The  response  spectra  for  strain  gage  2L  in  Shot  240-5, 
also  illustrated  in  Figure  16,  further  substantiates  that  the 
anticipated  displacements  were  not  observed  in  this  experiment. 
This  indicates  further  that  a  poor  detonable  gas  mix  probably 
existed  in  the  cavity  at  the  time  the  shot  was  fired. 
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The  displacement  response  spectra  obtained  from  the  dis¬ 
placement  history  reduced  from  strain  gage  2H  in  Shots  240-3c, 
240-4,  and  240-5  are  illustrated  in  Figure  17.  In  the  lower  fre¬ 
quency  range  where  these  response  spectra  are  most  valid,  the 
displacement  response  is  in  good  agreement  with  a  general  dis¬ 
placement  waveform  observed  in  these  experiments  and  the  dis¬ 
placement  response  spectra  obtained  from  the  low-frequency  ampli¬ 
fier  records  of  the  same  strain  gage.  Comparison  of  the  spectra 
obtained  for  Shot  240-4  with  that  from  240-3c  indicates  that 
their  amplitudes  in  the  10-  to  70-Hz  range  differ  only  by  a 
factor  of  approximately  2.3  as  opposed  to  the  factor  of  2.57 
observed  for  the  low-frequency  strain  gage  records.  This  des- 
crepancy  may  be  due  to  an  error  in  the  calibration  of  one  of  the 
strain  gage  amplifiers  used  in  the  experiments  or  to  some  other 
unknown  cause . 

The  displacement  response  spectra  obtained  from  the  reduced 
records  of  strain  gage  1H  in  Shots  240-6,  240-7  and  240-8  are 
shown  in  Figure  18.  Since  Shots  240-6  and  240-7  were  fired  at 
identical  initial  pressures  (220  psi )  ,  but  with  different  gas 
ratios,  these  spectra  should  provide  a  good  basis  for  evaluating 
the  importance  of  gas  mixtures. 

Shot  240-6  with  an  initial  methane-to-oxygen  ratio  of.lrl.o 
shows  a  significantly  lower  displacement  response  at  nearly  all 
frequencies.  In  the  10-  to  70-Hz  frequency  range,  the  displace¬ 
ment  response  for  Shot  240-7,  with  a  methane-to— oxygen  ratio  of 
1.2,  is  typically  1.5  times  greater  than  the  response  observed  in 
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Shot  240-6.  The  theoretical  ratio  in  energy  content  (per  unit 
volume)  fcr  a  1:1.5  mixture  compared  to  a  1.2  mixture  is  1:0.93 
(Reference  10).  The  data  of  strain  gage  1H  thus  appears  to  be 
quite  anomalous.  Careful  review  of  the  data  indicates  that  elec¬ 
trical  noise  during  Shot  240-6  made  the  strain  gage  baselines 
difficult  to  establish,  especially  for  the  close-in,  early 
recording  gage  locations.  Consequently,  an  error  in  baseline 
position  giving  smaller  late-time  (low-frequency)  displacements 
is  probably  responsible  for  the  difference  in  the  spectra  of 
strain  gage  1H  for  Shots  240-6  and  240-7.  However,  as  will  be 
discussed  later  with  reference  to  other  strain-gages,  displace¬ 
ment  spectra  for  these  two  tests  generally  show  larger  displace¬ 
ments  for  Shot  240-7  than  for  Shot  240-6. 

Shots  240-8  and  240-7  were  fired  with  the  same  methane/ 
oxygen  ratios  but  with  different  initial  pressures.  As  illus¬ 
trated  in  Figure  18  there  is  a  very  consistent  factor-of-two 
difference  between  the  displacement  response  for  these  two  shots 
for  all  scaled  frequencies  below  150  Hz.  Since  Shot  240-8  was 
fired  with  twice  the  initial  pressure  and  energy  as  Shot  240-7, 
the  doubling  of  displacements  indicates  a  fully  elastic  and, 
therefore,  completely  decoupled  behavior.  This  lack  of  increased 
coupling,  even  when  the  equilibrium  cavity  pressure  exceeds  the 
overburden  pressure  by  a  factor  of  two,  is  significant.  If  the 
salt  were  to  deform  plastically  as  was  observed  to  be  the  case 
for  the  natural  rock  salt  blocks  used  for  Shots  240-1  through 
240-5,  then  larger  displacements  would  have  definitely  been 
observed  at  the  higher  equilibrium  cavity  pressure  of  Shot  240-8. 
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Figure  16  Displacement  response  spectra  obtained  from  records 
of  strain  gage  2L  for  shots  240-3c,  240-4  and  240-5 
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Figure  19  Displacement  response  spectra  obtained  from 
records  of  strain  gage  2L  for  shots  240-6, 
240-7,  and  240-8. 
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Figure  20  Displacement  response  spectra  obtained  from  records 
of  strain  gage  2H  from  shots  240-6  to  240-9. 
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Figure  21  Displacement  response  spectra  obtained  from  large 
scale  calculations  illustrating  various  energy 
sources  and  rock  yield  models. 


if  the  salt  were  to  respond  purely  elastically  in  com¬ 
pression  but  were  to  fail  in  tension  at  very  low  stress  levels, 
then  equilibrium  cavity  pressures  up  to  three  times  the  over¬ 
burden  pressure  could  be  sustained  before  failure,  and  by  defi¬ 
nition,  increased  coupling  would  occur.  However,  the  well- 
established  failure  of  rock  under  conditions  of  high  pore  pres¬ 
sure  (Reference  11)  would  dictate  that  tensile  failure  at  lower 
cavity  pressures  should  probably  occur. 

The  artificially  pressed  salt  used  for  Shots  240-6  through 
240-9  might  be  expected  to  respond  elastically  to  compressional 
loading  due  to  the  work -hardened  nature  of  this  material.  It 
was  not  possible  to  anneal  the  salt  blocks  after  they  had  been 
pressed  and  the  high  pressures  (30,000  psi)  required  to  press 
the  salt  to  the  desired  density  indicated  that  severe  work¬ 
hardening  was  taking  place  during  the  pressing  operation.  The 
artificially  pressed  salt  still  had  some  porosity  and  permeabil¬ 
ity  and  it  is  surprising  that  pore-pressure  effects  at  the 
higher  firing  pressures  were  not  more  pronounced.  Additional 
efforts  to  investigate  this  phenomenon  would  certainly  be 
warranted. 


Excellent  strain  gage  displacement  records  were  also  obtained 
from  strain  gage  2  in  Shots  240-6,  240-7  and  240-8.  The  output 
of  this  gage  was  monitored  by  both  a  low-frequency  and  a  high- 
frequency  amplifier,  A  reasonably  good  high-frequency  record 
was  also  obtained  from  strain  gage  2  in  Shot  240-9.  The  displace¬ 
ment  response  spectra  calculated  from  these  strain  gage  records 
are  shown  in  Figure  19  for  the  low-frequency  amplifier  records 
in  Figure  20  for  the  high-frequency  amplifier  records. 
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The  displacement  response  spectra  shown  in  Figures  19  and 
20  generally  confirm  the  conclusions  drawn  from  the  response 
spectra  of  strain  gage  1H  shown  in  Figure  18.  Shots  240-8  and 
240-9,  fired  with  twice  the  initial  pressure  of  Shots  240-6  and 
240-7,  do  not  show  any  increased  coupling  due  to  the  higher 
equilibrium  cavity  pressures .  In  fact,  the  displacement  response 
spectra  for  strain  gage  2H  (Figure  20)  show  less  than  a  factor- 
of-two  increase  between  the  low-pressure  and  high-pressure  shots. 
This  less  than  expected  increase  in  displacement  response  can 
probably  be  attributed  to  the  poorer  recording  characteristics 
of  the  strain  gage  system  at  larger  distances  from  a  test  cavity. 

The  displacement  response  spectra  obtained  from  both  the 
low-  and  high-frequency  records  of  strain  gage  2  in  Shots  240-6 
and  240-7  do  not  show  the  large  discrepancy  which  appeared  in 
the  records  of  strain  gage  1H  as  already  discussed.  The  peak 
displacements  observed  in  Shot  240-7  are  still  slightly  higher 
than  the  240-6  displacements  although  theoretical  considerations 
(Reference  10)  indicate  otherwise.  Quite  possibly  the  assumption 
that  maximum  detonation  energy  release  is  realized  when  the 
reaction  products  are  just  CO  and  is  erroneous.  If  reactions 
leading  to  the  formation  of  CO2  can  contribute  to  the  energy  of 
detonation,  then  the  1:2  methane-oxygen  mixture  would  have  a 
higher-energy  output,  and  larger  displacements  would  be  expected. 

As  a  final  evaluation  the  results  of  the  spectral  analyses 
of  scaled  experimental  data  were  compared  with  displacement 
spectra  obtained  from  the  output  of  several  one-dimensional 
parametric  calculations.  These  calculations  were  recently  con¬ 
ducted  within  the  theoretical  portion  of  the  research  program 
(Contract  No.  DASA  01-70-C-0094)  that  covered  the  small-scale 
experiments  reported  upon  here.  These  calculations  were  designed 
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to  critically  demonstrate  and  evaluate  the  possible  effects  upon 
decoupling  phenomena  of  various  material  property,  explosive 
loading,  and  overburden  stress  conditions.  The  displacement 
response  spectra  for  several  calculations  particularly  comparable 
to  the  small-scale  experiments  are  illustrated  in  Figure  21.  A 
good  comparison  between  calculated  and  experimental  results  is 
shown  in  Figure  22.  Both  records  were  made  from  data  at  approxi¬ 
mately  4  cavity  radii  for  approximately  360  bar  cavity  equili¬ 
brium  pressure  and  180  bar  overburden.  Natural  cavity  resonance 
is  at  49  cps  (see  Section  4.2)  and  this  can  be  clearly  seen.  The 
calculated  amplitude  at  the  lower  frequencies  is  lower  than  the 
experimental  value  but  the  yield  surface  used  for  the  calculation 
may  be  higher  than  that  of  the  salt  used  in  the  experiment. 
Additionally  the  overburden  in  the  calculation  was  a  constant 
180  bars  and  in  reality  it  should  decrease  near  the  cavity.  Thus 
the  artificially  constant  overburden  in  the  calculation  sup¬ 
presses  yielding  near  the  cavity  and  results  in  the  lower  ampli¬ 
tudes  calculated  for  the  low  frequency  portion  of  the  section. 
Nevertheless  the  good  agreement  attests  strongly  to  the  validity 
of  the  experimental  technique,  the  strain  gage  diagnostic  system 
and  the  calculational  technique.  The  displacement  response 
spectra  obtained  from  the  data  of  strain  gage  2L  in  Shot  240-8 
(Figure  19)  is  also  very  comparable  to  two  of  the  calculated 
curves  shown  in  Figure  21.  Three  of  the  response  curves  are  for 
calculated  models  using  a  360-bar  equilibrium  cavity  pressure, 
equivalent  to  that  of  Shot  240-8.  One  of  these  curves  is  for  a 
model  with  a  180-bar  (2600  psi)  overburden  pressure — a  Mohr- 
Coulomb  yield  model  (MC-360  at  180).  The  response  spectra  for 
Shot  240-8  lies  intermediately  between  the  calculated  spectra  for 
the  Mohr- Coulomb  model  with  overburden  and  that  for  a  Von  Mises 
model  without  overburden  (VM-360  at  0).  This  comparison  indi¬ 
cates  that  some  yielding  of  the  salt  did  take  place  in  Shot  240-8. 
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Figure  22  Comparison  of  experimental  and  calculational  results 
for  180  bar  overburden,  approximately  360  bar  cavity 
equilibrium  pressure  at  4  cavity  radii. 


The  effects  of  having  a  very  weak  material  are  illustrated 
by  the  curve  for  a  Mohr-Coulomb  model  without  overburden  (MC-360 
at  0)  in  Figure  21.  Because  Mohr-Coulomb  yielding  is  pressure 
dependent,  the  absence  of  overburden  pressure  in  the  calculation 
allows  extensive  permanent  deformation  to  take  place.  This 
deformation  may  continue  into  late  times  as  radial  displacements 
effect  a  reduction  in  hoop  stresses  and  thus  mean  pressure,  which 
allows  for  additional  Mohr-Coulomb  yielding.  This  late-time 
yielding  is  illustrated  by  the  rapid  increase  in  displacement 
response  at  lower  frequencies.  Since  none  of  the  experimental 
results  show  such  low-frequency  behavior,  it  is  apparent  that  all 
of  the  salt  materials  used  in  the  experiments  were  not  subject  to 
such  late-time  yielding,  as  would  be  expected  with  the  presence 
of  overburden. 
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3.4  SUMMARY 


A  capability  for  conducting  small-scale  experiments  to  study 
decoupling  phenomena  was  successfully  developed,  tested,  and  used 
to  evaluate  some  of  the  basic  concepts  of  decoupling  theory. 
Because  of  the  low  cost  and  high  data  quality  associated  with 
small-scale  tests,  numerous  parameters  of  possible  interest  to 
decoupling  may  be  studied  in  a  cost  effective  fashion.  Because 
overburden  pressures,  test  materials,  and  explosive  loading  con¬ 
ditions  may  be  varied  as  desired  in  small-scale  experiments,  this 
capability  complements  both  field  tests  and  calculations  and 
would  be  preferable  to  in  situ  experiments  (both  large  and  small 
scale)  for  many  decoupling  studies. 

The  greatest  problems  that  had  to  be  overcome  in  developing 
the  small-scale  experimental  capability  were  associated  with 
detonating  me thane -oxygen  gas  mixtures  in  small  cavities  and  with 
eliminating  electrical  noise  from  the  strain  gage  diagnostic 
system.  Both  of  these  difficulties  were  overcome  with  the  suc¬ 
cessful  development  of  a  PETN  detonator  fired  by  a  long  delay 
line  of  MDF .  The  only  aspect  of  the  experimental  system  which 
could  benefit  from  further  development  is  the  strain  gage  diag¬ 
nostic  technique  for  measuring  radial  displacements.  Although 
high-auality  and  informative  displacement-time  histories  were 
obtained  with  the  strain  gage  technique,  it  is  felt  that  any 
future  studies  of  this  type  should  utilize  an  electromagnetic- 
velocity  gage  diagnostic  system. 
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Nine  small-scale  experiments  were  fired  using  the  small- 
scale  capability  developed  on  this  research  program.  With  the 
exception  of  electrical  noise  problems  in  two  of  the  experiments, 
excellent  strain  gage  data  were  obtained  from  gages  embedded  in 
the  salt  material  at  various  radii  from  a  2-inch"diameter  "de¬ 
coupling"  cavity.  Analysis  of  the  strain  gage  records  confirms 
the  basis  for  the  large-cavity  decoupling  theory.  In  natural 
rock  salt,  significant  increases  in  coupling  were  observed  as 
«_he  energy  released  by  detonation  of  a  methane-oxygen  mixture  in 
the  test  cavity  was  increased.  In  an  artificially  pressed  rock 
salt,  increased  coupling  was  not  observed  over  the  cavity  pres¬ 
sure  ranges  studied.  This  phenomenon  is  definitely  attributable 
to  the  work-hardened  nature  of  this  material  resulting  from  the 
cold  pressing  used  in  preparing  test  blocks. 

A  data  analysis  technique  involving  a  spectral  response 
analysis  was  developed  and  used  to  further  evaluate  the  data 
^ktained  from  the  small-scale  experiments .  The  displacement 
response  spectra  obtained  in  this  analysis  revealed  many  features 
of  the  data  that  were  not  apparent  in  the  reduced  displacement¬ 
time  curves  obtained  from  the  strain  gage  records.  The  displace¬ 
ment  response  spectra  confirmed  the  increased  coupling,  or  lack 
of  it,  observed  in  the  displacement-time  records.  Comparison  of 
these  response  spectra  with  those  obtained  from  calculated  (using 
a  one-dimensional  finite-difference  code)  displacement-time 
histories  revealed  very  good  agreement  between  the  two. 

Not  all  aspects  of  decoupling  phenomena  that  are  of  interest 
could  be  investigated  within  the  scope  of  the  program.  It  would 
have  been  particularly  informative  to  conduct  a  few  experiments 
(in  both  natural  and  artificially  pressed  rock  salt  test  blocks) 
using  a  higher  initial  gas  pressure  to  obtain  higher  explosive 
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energy  loadings  of  the  test  cavities.  The  results  of  such 
experiments  would  have  extended  the  range  and  degree  of  increased 
coupling  for  the  natural  rock  salt  tests  and  would  have  estab¬ 
lished  the  pressure  and  energy  levels  at  which  increased  coupling 
occurred  for  the  artificially  pressed  salt.  Any  future  small- 
scale  decoupling  experiments  should  probably  include  such  tests. 


SECTION  4 
CALCULATIONS 


4.1  CALCULATIONAL  TECHNIQUE 

4.1.1  Code  Capability.  Most  of  the  calculations  described 
in  this  section  were  carried  out  on  the  POD  code.  POD  is  a  one¬ 
dimensional  time-dependent  finite  difference  Lagrange  code,  which 
is  used  to  calculate  the  dynamic  motion  of  continue  including 
elastic-plastic  and  hydrodynamic  media.  The  differential  equa¬ 
tions  that  govern  the  motions  of  these  media  are  approximated  by 
difference  equations  applied  to  a  network  of  zones  that  describe 
the  physical  space  occupied  by  the  media.  Groups  of  contiguous 
zones  of  the  same  material  are  described  by  constitutive  equa¬ 
tions  which  are  coupled  to  the  equations  of  motions.  There  are 
three  possible  symmetries— plane ,  cylindrical  and  spherical.  The 
POD  code  offers  a  variety  of  initial  and  boundary  conditions,  an 
extensive  equation  of  state  and  yield  model  library  and  many 
other  specialized  options  developed  and  applied  on  a  wide  range 
of  problems  in  continuum  physics. 


The  two-dimensional  calculation  which  was  partially  com¬ 
pleted  in  the  course  of  this  contractual  effort  was  carried  out 
on  the  ELK  code  which  is  a  coupled  Eulerian-Lagrangian  code. 
This  code  was  developed  specifically  for  calculating  explosive 
phenomena  in  which  the  hydrodynamic  area  of  the  explosion  is 
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calculated  on  an  Eulerian  grid  and  the  elastic-plastic  regions 
away  from  the  explosion  are  calculated  on  a  Lagrange  grid 
coupled  to  the  Eulerian  grid.  There  are  tWo  possible  symmetries- - 
planar  and  axisymmetric .  Most  of  the  equation  of  state  and  yield  ■ 
models,  boundary  and  initial  condi tiohs,  and  other  special  options 
available  on  POD,  are  available  with  ELK. 

Both  POD  and  ELK  codes  are  compatible  with  Physics  Inter¬ 
national's  data  display  code  which  allow  time  histories,  profiles 
of  any  variable,  and  shock  trajectories  to  be  plotted  in  almost 
any  desired  format.  1 

i 

I 

In  the  calculations  discussed  in  this  report,  the  goal  was 
to  examine  the  stress  waves  generated  by  a  violent  explosion, 
after  they  had  propagated  for  seismologically  significant  dis¬ 
tances  into  the  surrounding  earth  media,  i  At  large  distance  the1 
kinetic  energy  density  of  outward  propagating  stress  waves 
quickly  decays  because  of  the  internal  friction  of  the  medium  and 
the  geometric  divergence  of  the  waves.  Beyond  a  characteristic 
"elastic  radius"  the  stress  waves  become  ordinary  acoustic  w^ves 
and  may  be  treated  by  the  theory  of  linear  elasticity.  It  is 
obvious  that  a  considerable1  saving  in  computer  !time  and1 storage  • 
space  can  be  realized  if  the  thermodynamic  iterations  and  plastic 
flow  computations  are  not  performed  for  those  finite-difference 
zones  that  lie  beyond  the  "elastic  radius"  of  the  energy  source 
of  the  problem. 

i  i 

Rather  than  modify  the  Lagrange  continua  mechanics  codes  to 

I 

skip  over  these  unnecessary  calculhtions  for  elastic  zones, 

.  .  ! 

Physics  International  developed  an  efficient  elastic  finite- 
difference  code  which  couples  to  the  Lagrange  codes.  This  code,  1 
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called  LEEK,  computes  elastic  zones  about  40  times  faster  than 
the  Lagrange  codes.  LEEK  is  a  sophisticated  elastic  code  and  may 
be  used  by  itself  as  a  tool  for  studying  seismic  wave  propaga¬ 
tion.  The  code  is  flexible  enough  to  compute  problems  in  which 
the  elastic  constraints  are  arbitrary  functions  of  position  and 
in  which  zones  are  arbitrary  quadrilaterals.  Either  stresses  or 
displacements  may  be  specified  as  the  boundary  conditions.  A 
description  of  the  LEEK  code  is  provided  in  Appendix  A. 

I 

I 

It  is  practical  with  the  aid  of  advanced  high  speed  com¬ 
puters  such. as  the  CDC  7600  to  use  the  LEEK  code  to  calculate 
detailed  problems  of  seismological  interest.  For  example,  a  one 
dimensional  LEEK  grid  with  1  meter  zone  width  and  10  ,000  zones 
(10  km)  can  be  coupled  to  an  elastic-plastic  Lagrange  grid  con¬ 
taining  an  explosive  or  fault  source.  Once  the  activity  of  the 
sjource  has  subsided,  the  Lagrange  grid  can  be  disconnected  and 
the  boundary  of  the  LEEK  code  can  be  maintained  at  the  equili¬ 
brium  stress  state,  while  the  stress  waves  generated  by  the 
source  are  observed  in  detail  as  they  propagate  out  several  kilo¬ 
meters.  A  two-dimensional  problem  utilizing  LEEK  would  be  of 
much  greater  interest  because  of  the  inclusion  of  free  surfaces 
and  layering.  However,  to  be  practical,  the  zonal  resolution 
.would  be  1  coarser ,  for  example  on  the  order  of  5  meters, and  the 
extent  of  the  grid  might  be  limited  to  three  or  four  kilometers. 
Stress,  displacement  and  velocity  histories  can  be  saved  at 
interesting  locations  and  used  as  gage  setting  information,  as 
data  prediction  or  interpretation  aids,  or  as  forcing  functions 
for  observing  structural  response.  The  salient  features  of  this 
type  of  calculation  are  the  inclusion  of  free  surfaces, 
arbitrary  layers  and  compatible  equation  of  state  models.  It 
should  be  noted  that  an  increase  in  zone-cycle  time  results  when 
a  simple  compaction  equation  of  state  is  used  rather  than  a 
linear  equation  of  state  in  the  LEEK  code. 
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4.1.2  Calculational  Procedure.  The  one-dimensional  calcu¬ 
lations  reported  herein  were  carried  out  in  two  basic  steps.  A 
hydrodynamic  source  calculation  was  made  and  the  pressure  history 
at  the  cavity  wall  was  saved  on  tape.  This  boundary  tape  was 
then  used  to  drive  a  Lagrange  grid  of  salt.  The  LEEK  grid  of 
purely  elastic  salt  was  coupled  to  the  outer  boundary  of  the 
Lagrange  grid  at  the  appropriate  time.  Stress  and  velocity 
histories  were  saved  on  tape  at  locations  of  1  cavity  radius  (the 
cavity  wall),  4  cavity  radii,  8  cavity  radii  and  10.7  cavity 
radii  (the  Lagrange — LEEK  interface).  The  cavity  radium  for  all 
the  ore-dimensional  calculations  was  1674  cm  corresponding  to  the 
estimated  volume  of  the  STERLING /MIRACLE  PLAY  cavitv.  The  LEEK 
grid  extended  from  10.7  cavity  radii  to  25.25  cavity  radii,  or 
nearly  half  a  km  from  the  source  center.  The  calculational  pro¬ 
cedure  is  illustrated  schematically  in  Figure  23. 

The  hydrodynamic  source  calculations  were  of  two  types--a 
mixture  of  methane-oxygen  detonated  at  the  center  of  the  cavity 
and  a  nuclear  source  (with  a  higher  shock  reverberation  frequency) 
also  detonated  at  the  center  of  the  cavity.  In  ail  the  source 
calculations  the  cavity  wall  was  fixed — the  very  small  effect  of 
the  wall  motion  on  the  gas  pressure  being  neglected. 

Three  detonable  source  calculations  (all  with  mixtures  of 

CH.  +  1.50-)  were  made,  corresponding  to  180  bars,  360  bars  and 
4  ^ 

720  bars  equilibrium  cavity  pressure.  All  three  calculations 
were  carried  out  to  0.1  sec. 

One  "nuclear"  calculation  was  made  as  representative  of  a 
high  frequency  source  with  the  same  equilibrium  cavity  pressure, 
360  bars,  as  one  of  the  gas  calculations.  This  high  frequency 
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Figure  23  Calculational  procedure  for  the  one-dimensional 
POD  problems. 


source  consisted  of  a  892  cm  diameter  ball  of  y  =  1.3  gas  it 

_  o  3 

1.21  x  10  gram/cm  and  18  kilobars .  The  rest  of  the  caviuy 
volume  was  void.  Since  the  void  occupied  about  98%  of  the  total 
volume,  the  temperature  of  the  expanded  fireball  was  high  and  the 
source  frequency  higher  (by  about  a  factor  of  3)  than  the  STERLING 
nuclear  source.  A  more  realistic  model  of  a  nuclear  source  such 
as  STERLING  would  require  the  cavity  volume  outside  of  the  fire¬ 
ball  to  be  filled  with  ambient  air  and  would  require  radiation 
diffusion  between  the  hot  fireball  gases  and  the  surrounding  cold 
air.  (A  summary  of  the  hydrodynamic  source  calculations  is 
given  in  Table  3.)  The  cavity  wall  pressure  histories  of  these 
source  calculations  were  then  applied  to  a  Lagrange  grid  of  salt 
(Figure  23) .  The  Lagrange  grid  was  chosen  large  enough  so  that 
the  stress  wave  reaching  the  outer  edge  was  always  elastic.  The 
grid  started  out  with  two-hundred  5  cm  wide  zones  and  was  rezoned 
in  stages  to  10  cm  wide  zones,  20  cm  wide  zones,  40  cm  wide  zones 
and  finally  80  cm  wide  zones.  Upon  the  fourth  rezone,  the  LEEK 
grid  of  three-hundred  80  cm  wide  zones  was  added  and  the  problems 
run  for  about  0.1  second. 

In  the  salt  calculations  several  parameters  were  varied. 

Three  different  gaseous  detonation  sources  and  one  high  frequency 
nuclear  source  were  used  to  drive  salt  with  two  different  yield 
models  (a  von  Mises  and  Mohr-Coulomb  model) ,  three  different 
equations  of  state  (one  linear  and  two  compaction) ,  and  over¬ 
burdens  of  zero,  180  bars  and  360  bars.  The  overburden  pressures 
were  generally  constant  throughout  the  Lagrange  grid,  however, 
one  calculation  was  carried  out  with  the  equilibrium  cavity 
stress  distribution  for  a  case  in  which  no  yielding  occurred 
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prior  to  the  explosion.  It  is  possible  to  read  in  equilibrium 
overburden  distributions  around  the  cavity  with  yielding  by 
utilizing  Physics  International's  Static  Code,  but  this  was  not 
carried  out  due  to  a  lack  of  time  and  money. 


A  summary  of  the  various  salt  calculations  is  given  in 
Tablf  4  and  a  description  of  the  equation  of  state  models  and 
yield  models  is  given  in  Appendix  B. 

Several  one-dimensional  calculations  were  carried  out  with 
various  zoning  schemes  and  stability  criteria  to  estimate  the 
effect  of  the  calculational  technique  on  the  answers  obtained. 
Problems  were  run  with  constant  80  cm  wide  zones  (rather  than 
rezoning  in  four  stages  from  5  cm  zones  to  80  cm  zones)  and  with 
constant  160  cm  wide  zones. 

These  same  problems  were  re-run  with  different  values  of  the 
safety  factor  for  the  Courant  stability  condition.  In  problems 
involving  the  propagation  of  elastic  waves  across  a  large  number 
of  finite  difference  zones,  artificial  dispersion  of  the  wave 
occurs.  This  artificial  dispersion  of  the  waveform  can  be 
minimized  by  maintaining  a  constant  zone  width  and  operating  with 
a  time  step  as  close  to  the  Courant  limit  as  possible,  without 
jeapordizing  the  stability  of  the  difference  equations.  This 
was  done  ii  all  the  calculations  discussed  above. 

In  the  two-dimensional  problem,  the  nuclear  source  and  salt 
grids  were  calculated  simultaneously.  The  cavity  for  this  calcu¬ 
lation  was  of  pillbox  geometry,  with  a  diameter  five  times  the 


SUMMARY  OF  ONE -DIMENSIONAL  SALT  CALCULATIONS 
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SUMMARY  OF  ONE-DIMENSIONAL  SALT  CALCULATIONS  (cont. 
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SUMMARY  OF  ONE -DIMENSIONAL 
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height  and  a  cavity  volume  equal  to  that  of  the  one-dimensional 
spherical  calculations.  The  source  was  a  984  cm  diameter  fire¬ 
ball  of  y  =  1.3  gas  at  1.21  x  10  ^  gram/cm^  and  15.6  kbars  sur¬ 
rounded  by  ambient  air.  However,  there  was  no  radiation  diffu¬ 
sion  mechanism  included  in  the  description  of  the  fireball  gas 
or  surrounding  air.  The  salt  was  described  by  a  linear  equation 
of  state  (P  =  0.265  y )  and  a  Mohr-Coulomb  yield  model  [Y  =  min 
(50  +  0.33  P,  200  bars)] ,  but  zero  overburden.  The  zoning  of  the 
Eulerian  source  grid  and  Lagrangian  salt  grid  is  shown  in 
Figure  58.  The  calculation  was  run  to  0.0075  sec,  a  time  when 
the  stress  wave  was  just  reaching  the  edge  of  the  Lagrange  grid. 
At  this  point,  when  the  LEEK  grid  was  to  be  added  on,  the  calcu¬ 
lation  was  terminated. 

4-1.3  Spectral  Analysis.  As  pointed  out  previously,  stress 

and  velocity  histories  were  recorded  at  1,  4,  8  and  10.7  cavity 

radii.  The  velocity  histories  were  subjected  to  a  spectral 

analysis  in  the  following  manner.  The  velocity  history  was 

integrated  to  yield  the  displacement  history.  The  displacement 

time  function  was  then  used  to  drive  an  undamped  mass-spring 

system  (Figure  45).  The  equation  of  motion  is  solved  and  the 

Displacement  Response  Spectra  (DSPEC)  and  Fourier  Velocity 

Amplitude  Spectra  are  computed.  From  the  Displacement  Response 

Spectrum,  the  Absolute  Acceleration  Response  Spectrum 
2 

[  ASPEC  =  w  DSPEC]  can  be  derived.  The  spectra  calculated  for 
the  one-dimensional  problems  are  valid  over  the  approximate  range 
10  cps  to  2,000  cps .  The  lower  limit  of  this  range  is  imposed  by 
the  problem  time  of  about  0.1  sec  and  the  upper  limit  is  governed 
by  the  onset  of  zonal  noise.  The  Fourier  Velocity  Spectra,  which 
is  always  equal  to  or  less  than  the  Displacement  Response 
Spectra,  is  not  strictly  accurate  as  the  driving  functions  are 
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terminated  before  they  die  away  naturally.  In  most  cases,  as 
will  be  discussed  below,  the  final  character  of  the  response 
spectra  is  established  by  the  time  the  driving  functions  are 
terminated.  A  detailed  discussion  of  the  various  response 
spectra  can  be  found  in  Reference  12. 

f 

4 . 2  RESULTS 

The  cavity  wall  pressure  histories  for  the  three  detonating 
gas  source  calculations  are  shown  in  Figures  24,  25  and  26.  The 
initial  gas  conditions  were  taken  from  Reference  13  and  an  ideal 
gamma  law  description  of  the  detonation  products  was  assumed. 

The  zoning  of  these  calculations  (80  zones  initially)  was  suffi¬ 
cient  to  resolve  the  Chapman-Jouguet  detonation  pressure  before 
the  wave  impinged  on  the  cavity  wall.  At  appropriate  intervals 
after  reflection  of  the  spherical  detonation  wave,  the  gas  cal¬ 
culation  was  dezoned  to  speed  up  the  calculation.  Representa¬ 
tive  pressure  profiles  for  the  three  source  calculations  are 
shown  in  'i.gures  27  and  28.  The  equilibrium  cavity  pressures 
for  the  three  calculations  are  180  bars,  360  bars  and  720  bars. 

The  cavity  wall  pressure  history  and  representative  pres¬ 
sure  profiles  for  the  nuclear  source  with  no  surrounding  air,  are 
shown  in  Figures  29  and  30.  The  equilibrium  cavity  pressure  is 
360  bars.  A  second  nuclear  source  calculation  with  ambient  air 
mixed  in  with  the  fireball  gas  was  made  and  the  wall  pressure 
history  is  shown  in  Figure  31.  The  fireball  and  ambient  gases 
were  artificially  pre-mixed  since  no  radiation  diffusion 
mechanism  is  presently  available  in  the  POD  code.  Artificially 
pre-mixing  the  gases  gives  the  proper  frequency  of  wall  pressure 
pulses  (Figure  31)  and  adjusting  the  initial  size  of  the  mixed 
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Figure  24  Cavity  wall  pressure  history  for  cavity  equilibrium 
pressure  of  180  bars  (CH.  +  1.509  +  1  atmos.  air  at 
300OK.  '  i  i 
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Figure  25 
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300  °K ) . 
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Figure  27  Selected  pressure  profiles  for  gas  source  calculation 
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Figure  28  Selected  pressure  profiles  for  the  gas  source  calcula¬ 
tions  . 
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Figure  29 
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Figure  30  Selected  pressure  profiles  for  the  nuclear  source 
calculation. 
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Figure  31  Cavity  wall  pressure  history  for  cavity  equilibrium 
pressure  of  360  bars  (nuclear  source  mixed  with  air 
in  cavity) . 


fireball  adjusts  the  size  of  the  initial  pressure  peak  measured 
at  the  wall.  However,  because  of  limitations  of  time  and  funds, 
the  second  nuclear  calculation  was  never  used  to  drive  a  salt 
calculation. 

The  frequency  of  pressure  pulses  in  the  detonating  gas  cal¬ 
culations  is  in  good  agreement  with  the  predicted  gas  dynamics  of 
the  MIRACLE  PLAY  series  (Reference  14)  ,  however  the  nuclear 
source  calculation,  with  no  air,  shows  a  considerably  higher 
frequency  of  pressure  pulses  than  shown  in  Reference  14.  The 
second  nuclear  source  calculation  with  mixed  fireball  and  ambient 
gases,  exhibits  a  frequency  of  wall  pressure  pulses  more  in 
agreement  with  those  predicted  for  the  STERLING  event  (Reference 

14)  . 

Representative  results  of  the  salt  calculations  are  shown  in 
Figures  32  through  44.  Radial  and  tangential  stress  histones  at 
13176  cm  or  8  cavity  radii  are  shown  in  Figure  32  for  the  case  of 
a  linear  equation  of  state,  180  bar  overburden,  Mohr-Coulomb 
yield  and  360  bar  cavity  equilibrium  pressure.  The  character  of 
these  histories  reflect  the  pulsed  character  of  the  source  pres¬ 
sure  history.  However,  the  character  of  the  salt  response 
changes  significantly  when  the  overburden  is  omitted  (for  the 
case  of  Mohr-Coulomb  yielding)  and  changed  even  more  so  when  a 
compaction  equation  of  state  is  used  (Figure  33).  As  wc uld  be 
expected,  the  overburden  suppresses  extensive  plastic  yielding 
around  the  cavity  and  the  hysteresis  associated  with  the  compac¬ 
tion  equation  of  state  modifies  the  wave  speeds  of  the  stress 
waves  and  relief  waves  tending  to  cause  successive  pulses  to 
coalesce.  The  effect  of  compaction  is  then  to  convert  energy 
from  the  higher  frequency  modes  to  lower  frequency  modes. 
Increasing  the  overburden  has  the  opposite  effect  of  shifting 
energy  to  higher  frequency  modes. 
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Figure  32 


Typical  calculated  stress  histories  in  the  SALT 
from  a  gaseous  detonation  in  a  spherical  cavity. 


Radial  stress  (bars) 


•  8  Cavity  radii  Time  (sac) 

•  Mohr-Coulomb  yield 

•  Equilibrium  cavity  pressure  =  360  bars 

Figure  33  Typical  radial  stress  histories  showing  the  effects 
of  overburden  and  compaction. 
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Figure  35  Velocity  histories  showing  the  effect  of  overburden 
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Figure  36  Velocity  histories  showing  the  effect  of  compaction 
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Figure  37  Velocity  histories  showing  the  effect  of  overburden 
and  yield  model. 
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Figure  38  Peak  particiQ  velocity  as  a  function  of  radius  show¬ 
ing  the  effect-  of  compaction  and  overburden. 
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Figure  40  Peak  particle  velocity  us  a  function  of  rac.ius 

showing  the  effect  of  cavity  equilibrium  pressure. 
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Figure  41  Peak  particle  velocity  as  a  function  of  radius  for 
nuclear  and  gas  explosions. 
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Figure  42  Displacement  histories  showing  the  effect  of  yield 
model  and  overburden. 
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Figure  44  Displacement  histories  for  gas  and  nuclear  source, 
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The  fundamental  cavity  frequency,  calculated  assuming  a 
cavity  in  an  infinite  elastic  media  is  given  by: 

*  _  I  (1  --2v)j.  %.  =  f§.  (c  } 

w  -  7  1  -  v  R  tiR  '  P 

where  v  =  Poisson's  ratio 

Cr  =  longitudinal  sound  speed 
Ij 

r  =  cavity  radius 

C  =  shear  wave  velocity 
s 

For  the  linear  equation  of  state  salt  used  in  the  calculations: 

v  =  0.2713 
C.  =  0.4bl0  cm/psec 

Ij 

r  =  1647  cm 
w*  =  48.8  cps 

The  cavity  "ringing'  is  highly  damped  because  of  elastic  energy 
radiating  away  and  spherical  divergence.  However,  the  ringing 
can  be  seen  clearly  in  the  stress  history  in  the  salt  driven  by 
the  high  frequency  nuclear  source  (Figure  34).  Although  a  Mohr- 
Coulomb  yield  model  was  used,  there  was  sufficient  overburden  to 
suppress  most  of  the  plastic  yielding  so  that  the  cavity  ringing 
can  clearly  be  seen.  It  is  less  obvious  in  the  equivalent 
gas  calculation  primarily  because  the  frequency  of  pressure  pulses 
(62  cps)  is  very  close  to  the  natural  cavity  frequency.  However, 
as  will  be  described  later  the  fundamental  cavity  frequency  can  be 
seen  unmistakeably  in  the  response  spectra. 

Selected  particle  velocity  histories  also  at  8  cavity  radii 
(13176  cm)  are  shown  in  Figures  35,  36  and  37.  Figure  35  portrays 
the  effect  of  overburden  with  a  Mohr-Coulomb  yield  model  and 
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Figure  36  shows  the  effect  of  varying  the  degree  of  compaction. 

A  comparison  of  yield  models  is  presented  in  Figure  37,  where  it 
is  seen  that  the  particle  velocity  for  the  case  of  Von-Mises 
yield  with  no  overburden  is  very  similar  to  that  of  Mohr-Coulomb 
yield  with  overburden.  Yielding  occurs  at  much  lower  stresses 
with  the  Mohr-Coulomb  model  and  is  suppressed  when  overburden  is 
included.  Particle  velocity  histories  for  the  other  one- 
dimensional  salt  calculations  are  given  in  Appendix  D. 


Plots  of  peak  particle  velocity  as  a  function  of  radius  are 
given  in  Figures  39  through  41.  As  expected  the  peak  particle 
velocity  is  inversely  proportional  to  radius  when  the  wave  is 
elastic,  in  all  cases  the  extent  of  plastic  yielding  was  well 
within  the  LEEK  interface  at  17,647  cm  (10.7  cavity  radii).  The 
variation  of  peak  velocity  near  the  source  is  more  complicated 
because  of  yielding  and  compaction.  A  close  examination  of 
these  plots  shows  the  filtering  effects  of  compaction  and 
yielding.  For  example,  a  compactible  equation  of  state  results 
m  reduction  in  peak  velocity  by  a  factor  of  approximately  2.3 
(Figure  38)  while  extensive  Mohr-Coulomb  yielding  reduces  peak 
velocity  by  a  factor  of  about  1.4  (Figure  39)  over  that  of 
simple  Von  Mises  yielding.  Doubling  the  cavity  equilibrium 
pressure  with  and  without  overburden  increases  the  peak  velocity 
by  a  factor  of  approximately  1.6  (Figure  40).  As  will  be  dis¬ 
cussed  later,  the  maximum  displacement  and  low  frequency  response 
spectrum  is  increased  and  the  peak  particle  velocity  is  decreased 
when  the  media  yields  and/or  is  compacted. 


Figure  41  compares  the  peak  particle  velocity  versus  radius 
behavior  for  a  gas  and  nuclear  source.  They  are  remarkably 
similar  except  for  an  anomalous  dip  in  the  peak  velocity  at 
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about  10  cavity  radii  from  the  nuclear  source.  It  is  suspected 
that  this  is  a  result  of  a  loss  in  resolution  of  the  tip  of  the 
very  narrow  spike  of  the  stress  wave.  This  anomalous  dip  would 
not  have  occurred  if  the  zoning  in  the  salt  were  finer. 

Typical  displacement  histories  at  10.7  cavity  radii 
(17.647  cm)  are  shown  in  Figures  42,  43  and  44.  These  show  the 
increased  decoupling  with  increased  overburden  (Figure  42)  and 
demonstrates  the  importance  of  including  overburden  in  calcula¬ 
tions  or  experiments  of  decoupling  phenomena.  In  the  case  of  no 
overburden  and/or  compaction  (Figure  43) ,  the  displacement  con¬ 
tinues  to  increase  at  the  termination  of  the  calculation  and 
suggests  that  such  calculations  if  carried  out  to  much  longer 
times  would  demonstrate  significant  seismological  effects  of 
yield  model  and  compaction  model  variation. 

The  comparison  of  displacement  histories  generated  by  gas 
and  nuclear  source  (Figure  44),  shows  that  the  nuclear  source, 
which  was  fully  decoupled  (almost  no  yielding  near  the  cavity  was 
observed)  results  in  peak  displacement  nearly  an  order  of  magni¬ 
tude  less  than  that  resulting  from  the  gas  source.  Part  of  this 
may  be  due  to  a  tendency  of  the  cavity  to  resonate  with  the 
gaseous  detonation  pressure  pulses  as  the  two  frequencies  are 
very  close  (49  cps  resonant  frequency  excited  by  62  cps  pressure 
pulses).  This  comparison  demonstrates  the  importance  of  source 
frequency  on  decoupling  behavior. 

Several  calculations  were  made  to  investigate  the  effects  of 
modifying  the  zoning  and  Courant  stability  condition  safety 
factor.  With  coarser  zoning  the  wave  fronts  were  spread  out  but 
not  seriously.  Coarser  zoning  would  be  quite  adequate  for 
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extended  calculations  where  seismologically  significant  effects 
are  being  studied.  Even  the  higher  frequency  components  of  the 
waves  (several  hundred  cps)  do  not  appear  to  be  seriously 
affected  by  coarser  zoning  as  will  be  shown  in  the  spectral 
analysis  given  later.  With  a  high  Courant  stability  condition 
safety  factor  there  is  significant  noise  generated  in  the  relief 
waves.  This  is  caused  by  an  increased  interpolation  error  in  the 
difference  equations  as  a  result  of  the  larger  time  step.  It  is 
therefore  felt  that  the  advantage  of  reduced  numerical  dispersion 
of  the  stress  wave  gained  by  using  a  very  high  safety  factor  is 
outweighed  by  the  excessive  noise  generated  in  the  relief  wave. 

4.3  DATA  ANALYSIS 


The  velocity  histories  monitored  for  all  the  calculations 
were  subjected  to  a  spectral  analysis.  The  velocity  histories 
were  integrated  to  give  displacement  histories  whereupon  the  dis¬ 
placement  response  spectra  and  Fourier  amplitude  spectra  were 
calculated  (see  Figure  45  and  Reference  12  for  example) . 

Examples  of  displacement  response  spectra  are  shown  in  Figure  46 
where  the  highlights  such  as  theoretical  fundamental  cavity  fre¬ 
quency,  source  shock  reverberation  frequency  and  the  spectral 
window  (determined  by  pulse  length  and  onset  of  zonal  noise)  are 
illustrated.  Zonal  noise  appears  to  be  noticeable  in  the  spectra 
at  2,000  to  2,500  cps.  The  transit  time  for  a  sound  wave  across 
a  single  zone  is  80  cm/0.4510  ysec  =  177  ysec,  which  corresponds 
to  5,650  cps,  so  the  observed  onset  of  noise  appears  to  be 
related  to  2  to  3  zone  transit  times.  Approximately  8  higher 
modes  of  the  fundamental  shock  reverberation  frequency  are 
resolvable  before  the  onset  of  noise. 


Ill 


Figure  45  Mass-spring  system  used  to  determine  the  Fourier 
and  response  spectra. 
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Figure  46  Structure  of  the  displacement  response  spectra. 


113 


i 


I 


I  i  !  I  |  I  J 

' 

Figure  47  compares  the  response  spectra  of  five  1/2  pulses  with 
that  of  the  first  pulse  only  for  a  case  with  little  yielding  and 
no  compaction.  It  is  seen  that  the  first  pulse  establishes  the 
basic  frame  of  the  response  spectra  and  determines  the|  maximum 
displacement.  Successive  pulses  build  up  the  primary  and  higher' 
order  maxima  end  also  excite  some  ,natural  cavity  resonance. 

J  .  i 

1 

I 

The  salt  response  for  thh  high  frequency  nuclear  source  is 

,  *  j 

compared  to  that  of  the  corresponding  gas  source  in  Figure  48  and 

«  > 

it  is  seen  that  the  spectrum,,  especially  the  low  frequency 
portion,  of  the  nuclear  source  is  lower  in  amplitude  (showing 
better  decoupling).  The 1  pressure  ,pulse  frequency  for  the  nuclear 
source  is  about  430  cps ,  while  that  of  the  gas  source  is  62  cjbs. 

The  excitation  of  natural  cavity  resonknce  is  clearly  evident 
for  !  the  nuclear  source  and  is  readily  apparent  for  the!  gas 

'  !  1  .  i 

source . 

:  .  l 

1 

The  effect  of  radius  on  the  response  spectra  is  shown  in  ,  i 
Figure  49  for  a  case  with  little  yielding  and  no  compaction  and 
all  but  the  higher  frequency  content  seems  to  scale  as  expected; 

i  1 

The  effect  of  varying  overburden  with  a  Mohr-Coulomb  yield  i 
model  is  shown  in  the  response  spectra  of  Figure  50.  Increasing 
the  overburden  clearly  suppresses  the  peak  displacement  (uj  =  o,  . 

The  amplitude  of  the  fundamental  shock  reverberation  frequency 
(~  62  cpc)  increases  with  overburden.  As  the  overburden  is  i 

i  i  •  ’  1  .  _  i 

removed  more  plastic  yielding  can  pccur  and  energy  shifts  from 
higher  frequency  to  lower  frequency  modes. 

1  i  , 

i  i  ; 
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Figure  47  Response  spectra  for  the  first  pulse  and  for  the  first 
five  and  one  half  pulses  calculated  for  an  equilibrium 
cavity  pressure  of  360  bars,  overburden  of  180  bars, 
Mohr  Coulomb  yield  at  4  cavity  radii. 
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Figure  48  Displacement  response  spectra  showing  the  effect  of 
source  frequency. 
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Figure  49  Displacement  response  spectra  showing  the  effect  of 
distance  from  the  source. 
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Figure  50  “^placement  response  spectre  showing  the  effect 


The  effects  of  varying  overburden  with  a  von  Mises  yield 

burden ^h  ^  “  ,trlMn9  <,i’Uro  “>  3in~  "ith  -  without  over- 
burden  there  was  not  extensive  yielding  in  either  of  the  two 
cases  calculated. 


with  '  rmPariS°n  °f  0,0  V°”  MiS6S  and  M°hf Coulomb  yielding 
With  overburden  (Figure  51)  is  interesting.  All  things  being 

qual,  it  would  be  expected  that  the  amplitude  of  t-e  low  fr  - 
quency  spectrum  would  be  higher  with  the  Mohr-Coulomb  model 
since  more  yielding  occurs.  From  Figure  51,  the  opposite 
Wears  to  be  true.  The  reason  is  that  in  the  «ohr-Coulomb  cal¬ 
culation  the  overburden  was  constant  180  bars,  starting  from  the 

ZTol  mw  T  OXt°ndlng  OUtWard‘  ln  thC  calculation 

e  equilibrium  overburden  distribution  was  put  in,  starting  at 

0  bars  at  the  cavity  wall  and  increasing  to  180  bars  at  very 

arge  distances.  The  stress  reached  170  bars  at  about  2-1/2 

cavity  radii.  Evidently  enough  yielding  occurred  close  to  the 

cavity  to  cause  the  peat  displacement  to  be  higher  than  in  the 

tohr-Cou  omb  calculation.  This  demonstrates  that  including  the 

equilibrium  stress  distribution  around  the  cavity  does  have  a 

detectable  effect  and  should  be  included  for  accurate  calculation 
of  decoupling  behavior. 

Figure  52  portrays  the  response  spectra  showing  the  effects 
of  compaction,  with  compaction  there  is  a  dramatic  shift  of 
energy  from  high  to  low  frequency  modes  and  the  shock  reverbera¬ 
tion  frequency  maxima  are  obscured  or  completely  missing,  if 
overburden  were  included  in  the  compaction  calculations,  the 
maximum  displacements  (w  -  o)  observed  would  be  lower,  however 
the  spectra  would  still  not  show  the  reverberation  frequency  ' 
maxima  and  there  would  still  be  a  substantial  shift  of  energy 
from  higher  frequency  to  lower  frequency  modes. 


*.19 


Displacement  response  spec 


0.01 


0.001 


0.0001 


Figure 


von  Mises  yield 
overburden  =  0 

von  Mises  yield 
overburden  =  180  bars 

Mohr-Coulomb  yield 
overburden  =  180  bars 


10.7  Cavity  rauii 
Cavity  equilibrium  pressure 
®  360  bars 
P  -  0.265u 


.0 


-J _ I _ 

10.0  100.0 

Frequency  (cps) 


_ L_ 

1000.0 


51  Displacement  response  spectra  showing  the  effect 
of  overburden  and  yield  model. 
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Figure  52  Displacement  response  spectra  showing  the  effects 
of  compaction  and  yield  model. 


Figure  53  shows  the  effect  of  increasing  cavity  equilibrium 
pressure  on  the  displacement  response  spectra  with  and  without 
overburden.  With  the  higher  cavity  equilibrium  pressures,  more 
yielding  occurs  and  there  is  less  decoupling  as  expected. 

A  summary  of  maximum  displacements  for  most  of  the  calcula¬ 
tions  is  shown  in  Figure  54  as  a  function  of  radius.  The 
reduced  decoupling  effectiveness  of  yielding  and/or  ccmpactible 
media  is  evident  in  this  plot.  Compaction  and  lacV-  of  overburden 
with  Mohr-Coulomb  yielding  resulu  in  strongly  coupled  cavity 
response,  and  the  high  frequency  nuclear  source  is  seen  to  be  the 
best  decoupled  case. 

The  response  spectra  for  various  zoning  and  Courant 
stability  condition  safety  factors  are  shov/n  in  Figure  55.  Some 
losses  are  associated  with  coarser  zoning  but  the  effects  are 
minor.  Coarser  zoning  than  used  in  the  present  calculations 
would  not  seriously  affect  the  results  of  the  calculations  and 
would  save  computer  time,  especially  for  extended  runs. 

For  the  case  of  Mohr-Coulomb  yield,  overburden  and  linear 
equation  of  state,  the  response  and  Fourier  spectra  are  compared 
in  Figure  56 .  As  expected  the  Fourier  spectrum  is  equal  to  or 
less  than  the  response  spectrum.  The  spectra  coincide  at  the 
shock  reverberation  frequency  and  its  higher  modes  and  nearly 
coincide  at  the  natural  cavity  resonant  frequency.  The  maximum 
displacement  (which  from  Figure  47  results  from  the  first  pulse) 
is  larger  than  the  Fourier  amplitude  (which  is  the  maximum  dis¬ 
placement  after  the  entire  signal). 
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Displacement  response  spectra  (cm) 


Figure  53 


Displacement  response  spectra  showing  the  effects  o 
equilibrium  cavity  pressure  and  overburden. 
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Figure  55  Displacement  response  spectra  showing  the  effect  of 
zoning  and  the  courant  stability  condition  safety- 
factor. 


In  the  case  of  a  compactib]e  media  (Figure  57)  ,  the  response 
and  Fourier  spectra  virtually  coincide  since  the  displacements 
were  still  increasing  (Figure  43)  when  the  calculation  was 
terminated. 

For  an  explosion  source  (the  first  pressure  pulse  is  the 
largest) ,  the  degree  to  which  the  response  and  Fourier  spectra 
coincide  is  a  measure  of  the  compaction  and  yielding  or  other 
non-linear  behavior  that  occurs  in  the  media. 

4.4  PARTIAL  RESULTS  OF  THE  2-D  CALCULATION 

The  zoning  of  the  2-D  Lagrange  grid  is  shown  in  Figure  58. 
The  cavity  is  a  pillbox  with  a  diameter  five  times  its  height 
and  the  interior  of  the  cavity  is  represented  by  an  Eulerian 
grid.  The  initial  configuration  of  the  fireball  is  shown  in 
Figure  58.  The  fireball  is  approximately  a  984  cm  diameter 
sphere  of  y  =  1.3  gas  at  1.21  gram/ cm3  and  15.6  kbars  and  is 
surrounded  by  air  at  atmospheric  conditions. 

As  the  fireball  expands  it  first  strikes  the  cavity  directly 
above  it  and  spreads  rapidly  to  fill  the  entire  cavity.  The  gas 
dynamics  are  very  complicated  but  because  of  the  coarse  zoning  of 
the  Eulerian  grid,  most  of  the  finer  detail  was  lost.  Velocity 
histories  in  the  surrounding  salt  along  the  X  and  Y  axes  (Figures 
60  and  61)  show  that  at  a  given  distance  from  the  origin  the 
stress  wave  arrives  first  along  the  Y  axis  (the  fireball  expands 
at  a  higher  velocity  than  stress  waves  move  in  the  media)  and 
that  the  initial  wave  is  stronger  along  the  Y  axis. 
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Figure  57  Comparison  of  response  and  Fourier  spectra  for 
compactible  media  with  yielding  and  gas  source 
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Figure  58  Lagrange  grid  at  t  =  7500  usee.  Also  shown  are  locations  of 

stress  and  velocity  histories  and  initial  fireball  configuration. 


The  velocity  vector  plot  at  t  =  0.0075  sec  (Figure  59) 
indicates  that  the  pillbox  geometry  results  in  substantial  shear 
wave  energy  although  it  is  too  early  tc  tell  if  it  is  seismo- 
logically  important.  This  typ<i  of  plot  is  especially  useful  in 
visualizing  the  waves  as  they  radiate  into  the  salt  and  shows  the 
shear  waves  quite  clearly. 

The  X,  Y  and  shear  stress  histories  at  45°  from  the  axis  of 
symmetry  are  shown  in  Figures  62,  63  and  64  and  the  peak  X  stress 
is  seen  to  be  about  1/3  larger  than  the  peak  Y  stress. 

The  calculation  was  terminated  just  before  the  LEEK  grid  was 
to  be  added.  It  was  intended  to  run  the  calculation  to  t  =  0.05 
and  compare  the  stress  and  velocity  histories  with  those  of  the 
one-dimensional  calculations. 

4.5  SUMMARY 

The  calculations  presented  in  this  report  demonstrate  a 
technique  whereby  finite  difference  codes  may  be  used  to  study 
decoupling  phenomena.  Hydrodynamic  and  elastic-plastic  grids  are 
used  to  calculate  the  details  of  the  explosion  source  and  the 
complex  response  of  the  media  in  the  vicinity  of  the  cavity.  A 
purely  elastic  code,  LEEK,  is  coupled  to  the  elastic-plastic 
Lagrange  grid  to  permit  stress  waves  to  propagate  far  out  into 
the  media.  It  is  also  feasible  to  carry  out  near  teleseismic 
calculations  (~  3-5  km  from  the  source)  by  using  two-dimensional 
codes  coupled  to  a  LEEK  grid.  Such  calculations  wo'  Id  include 
die  effect  of  layering,  free  surfaces  and  material  compaction 
and  would  have  zonal  resolution  on  the  order  of  5-meters.  When 
run  on  advanced  computer  systems  such  as  the  CDC  7600,  the  cost 
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Figure  59  Velocity  vector  plot  of  the  _.agrange  grid  at  t  =  7500  usee 
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Figure  60  Velocity  (X)  history  parallel  to  the  axis  of 

symmetry  of  the  cavity  at  X  =  2412  cm  and  Y  =  0  err 
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Figure  63  Stress  (tyy)  history  perpendicular  to  the  axis  of 

symmetry  of  the  cavity  at  X  =  1998  cm  and  Y  =  1839  cm 
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of  these  calculations  would  be  reasonable.  The  results  of  this 
type  of  calculation  would  provide  insights  for  interpreting  near 
teleseismic  data  gathered  on  many  recent  underground  tests. 

In  the  course  of  the  one-dimensional  calculations  reported 
herein,  the  effects  of  yield  model,  overburden,  compaction, 
source  amplitude  and  frequency  were  studied.  An  evaluation  of 
the  salt  response  in  the  time  and  frequency  domain  revealed  that 
several  of  the  above  parameters  could  have  a  controlling  effect 
on  tne  seismic  signal  generated  by  the  explosion. 

For  a  given  cavity  and  explosion  source,  the  plastic 
yielding  in  the  vicinity  of  the  cavity  (in  the  first  few  cavity 
radii),  shifts  energy  from  the  high  frequency  modes  (tens  and 
hundreds  of  cps)  to  the  low  frequency  modes  which  are  able  to  be 
transmitted  to  teleseismic  distances.  Although  there  is  a  net 
loss  of  enerc/y  associated  with  yielding,  there  is  an  amplifica¬ 
tion  of  energy  in  the  low  frequency  modes . 

The  effect  of  increasing  overburden  is  to  suppress  yielding 
and  thereby  suppress  the  transfer  of  energy  to  the  lower  fre¬ 
quency  modes.  Equilibrium  overburden  stress  distribution  around 
the  cavity  is  important  in  the  calculation,  since  there  is  sig- 
yielding  in  the  immediate  area  of  the  cavity  which  has 
an  important  role  in  determining  the  .Low  frequency  content  of  the 
energy  spectrum. 

One  of  the  most  dramatic  effects  observed  in  the  calculations 
was  that  of  material  compaction.  Compaction  causes  a  large 
amount  of  net  energy  loss,  but  there  is  still  a  net  amplification 
in  the  low  frequency  content  of  the  energy  spectrum.  No  calcula¬ 
tions  of  a  highly  coupled  source  in  a  yielding,  compactible  media 


137 


with  overburden  was  made.  However,  it  is  not  clear  that  in  such 
a  case  that  the  energy  losses  in  the  yielding,  compacting  media 
would  be  large  enough  to  overcome  the  amount  of  energy  shifted 
from  the  higher  frequency  modes  to  the  lower  (and  seismically 
important)  modes.  The  object  of  decoupling  is  to  minimize  the 
low  frequency  amplitude  of  the  signal  and  the  calculational 
technique  used  in  the  present  study  can  be  a  valuable  tool  in 
studying  presently  held  views  of  large  cavity  decoupling  or  in 
optimizing  source,  cavity  and  media  properties  to  minimize  the 
low  frequency  signal  amplitudes  in  future  underground  tes^s. 

Most  of  the  calculations  reported  herein  were  made  for  a 
detonating  gas  source  relevant  to  the  MIRACLE  PLAY  Series.  One 
high  frequency  source  representative  of  a  nuclear  explosion  in  an 
evacuated  cavity  is  presented.  The  effect  of  increased  source 
frequency,  all  other  things  being  equal,  was  a  marked  reduction 
in  coupling  to  the  surrounding  media.  The  source  frequency  of 
the  nuclear  calculation  was  well  above  the  natural  cavity 
resonance  frequency  and  there  was  very  little  cavity  yielding . 
Consequently,  most  of  the  radiated  energy  resides  in  the  high 
frequency  portion  of  the  spectrum  (several  hundred  cps)  and  will 
decay  quickly. 

A  two-dimensional  calculation  to  investigate  the  effect  of 
cavity  geometry  on  the  observed  signal,  was  partially  completed 
and  indicated  that  substantial  shear  wave  energy  is  generated 
although  the  calculation  was  not  carried  out  far  enough  to 
determine  if  this  was  of  seismic  importance. 
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SECTION  5 
CONCLUSIONS 


A  capability  for  conducting  small  scale  experiments  with 
various  materials,  overburdens  and  sources  has  been  successfully 
developed.  Used  in  conjunction  with  field  tests  and  calcula¬ 
tions,  this  technique  can  be  used  to  isolate  and  study  parameters 
affecting  decoupling  and  can  lead  to  in  leased  understanding  of 
presently  available  data,  and  possibly  reduce  the  number  of 
underground  tests  required  in  the  future.  The  small  scale 
experiments  reported  herein  concenJ  rated  on  methane-oxygen 
detonations  in  spherical  cavities  in  salt.  The  salt  behavior  was 
in  general  agreement  with  large  cavity  decoupling  theory. 
Increased  coupling  was  observed  at  higher  pressures  in  natural 
salt  but  was  not  observed  over  the  pressure  range  studied  in 
artificially  pressed  salt  where  work-hardening  and  overburden 
apparently  suppressed  extensive  yielding. 

High  quality  data  was  obtained  using  a  strain  gauge  diagnos¬ 
tic  system,  however  even  better  data  can  oe  obtained  with  elec¬ 
tromagnetic  velocity  gage  systems  now  available  at  Physics 
International.  Subjecting  the  data  to  a  response  spectra 
analysis  further  clarifies  the  decoupling  effects  of  yielding  and 
overburden  as  it  provides  a  measure  of  the  transfer  of  energy 
from  higher  to  lower  frequency  modes.  The  salt  response  observed 
in  these  tests  is  also  in  general  agreement  with  the  response 
predicted  by  a  series  of  one-dimensional  finite  difference  cal¬ 
culations  . 
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The  experimental  apparatus  now  exists  to  perform  similar 
experiments  with  different  media  such  as  granite  or  tuff  and  is 
suitable  for  any  chemical  explosive  souice.  In  addition  to 
simulating  overburden,  the  apparatus  car  be  used  to  investigate 
non-spherical  cavities  and  evacuated  cavities  thereby  allowing  a 
study  of  varying  source  frequency.  With  an  evacuated  or  helium 
filled  cavity  it  is  possible  to  provide  a  source  with  amplitude 

and  frequency  characteristics  approaching  those  of  a  nuclear 
source . 

A  very  interesting  test  would  be  to  compare  a  given  source 
in  a  large  decoupled  cavity  in  an  elastic  media  with  the  same 
source  in  a  highly  coupled  cavity  in  a  yielding  compactible 
media.  A  response  spectra  analysis  of  the  velocity  history  at 
6  cavity  radii  would  permit  a  direct  measurement  of  the  low 
frequency  content  of  the  energy  spectrum  and  would  help  determine 
which  approach  might  be  better  for  effective  decoupling. 

A  calculational  technique  has  been  successfully  developed 
which  has  the  potential  of  solving  problems  in  the  near  tele- 
seismic  range  (several  kilometers).  This  technique  utilizes 
traditional  finite  difference  codes  to  calculate  the  explosion 
and  complex  response  of  the  media  near  the  cavity  and  couples 
this  to  the  purely  elastic  LEEK  code,  which  extends  the  calcula¬ 
tion  to  the  near  teleseismic  regime.  This  method,  which  was 
demonstrated  in  a  series  of  one-dimensional  calculations,  is  also 
applicable  to  existing  two-dimensional  codes,  and  when  used  on 
advanced  computer  systems  can  calculate  at  a  reasonable  cost 
problems  of  current  interest  such  as  the  near  teleseismic  data 
gathered  on  the  Diamond  Mine  and  Diamond  Dust  experiments . 
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In  the  one-dimensional  calculations  discussed  in  this 
report,  it  is  pointed  out  that  the  details  of  the  waves  propa¬ 
gating  in  the  overburdened,  yielding  and  compacting  media  near 
the  cavity  are  critically  important  in  determining  the  seismic 
content  of  the  signal.  For  example  in  a  highly  yielding  and 
compacting  media  there  can  be  a  large  energy  loss  in  the  media 
but  still  a  net  amplification  of  energy  in  the  low  frequency 
portion  of  the  spectrum  which  is  of  great  importance  in 

decoupling. 

Both  the  experimental  and  calculational  capabilities 
reported  herein  are  now  developed  to  the  point  where  they  can 
be  used  immediately  to  study  problems  of  specific  interest. 
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APPENDIX  A 
THE  LEEK  CODE 


In  a  typical  problem  of  cratering  of  a  solid  material, 

violent  hydi  jdynamic  motion  and  plastic  flow  is  confined  to  a 
relatively  i-mall  region  around  the  explosion  or  impact.  The 
kinetic  energy  density  of  outward  propagating  stress  waves 
quickly  decays  because  of  the  internal  friction  of  the  medium 
and  the  geometric  divergence  of  the  waves.  Beyond  a  charac¬ 
teristic  "elastic  radius,"  Re ,  the  stress  waves,  which  may 
still  be  of  kilobar  magnitudes,  become  ordinary  acoustic  waves 
and  may  be  treated  by  the  theory  of  linear  elasticity.  The 

elastic  radius  in  a  typical  cratering  problem  might  be  a  few 
crater  radii,  for  example.  In  a  two-dimensional  calculation 
of  cratering,  a  large  fraction  of  the  computational  time  in 
each  time  step  is  devoted  to  performing  thermodynamic  and 
plastic  flow  computations.  It  is  obvious  that  a  considerable 
saving  in  computer  time  and  storage  space  can  be  realized  if 
the  thermodynamic  iterations  and  plastic  flow  computations  are 
not  performed  for  those  finite-difference  zones  that  lie  beyond 
the  elastic  radius  of  the  energy  source  of  the  problem. 

Rather  than  modify  our  two-dimensional  Lagrange  continuum 
mechanics  code  to  skip  over  these  unnecessary  calculations  ior 
elastic  zones,  Physics  International  has  developed  an  efficient 
and  flexible  elastic  finite-difference  code  which  couples  to 
our  Lagrange  code.  This  code  is  called  LEEK  (Linear  Elastic 
Escape  KondiLion)?  it  computes  elastic  zones  about  40  times 

faster  than  our  Lagrange  code. 
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Because  of  its  computational  efficiency,  LEEK  greatly 
broadens  the  scope  of  problems  which  our  continuum  mechanics 
codes  can  economically  and  accurately  compute.  As  an  example, 
in  the  time  it  takes  for  a  crater  to  form,  stress  waves  may 
travel  on  the  order  of  10  crater  radii  from  the  point  of  impact. 

To  compute  the  full  crater  formation,  one  must  have  finite- 
difference  zones  out  to  about  six  crater  radii,  so  that  wave 
reflections  from  the  outer  boundary  will  not  return  to  the 
region  of  the  crater  during  the  time  of  its  formation  (Figure  A-l) . 
Such  reflections  would  be  unrealistic  anu,  duo  to  their  geometric 
convergence ,  could  ruin  the  computation.  It  would  be  prohib¬ 
itively  expensive  to  use  the  Lagrange  code  to  compute  all  of 
the  zones  of  this  problem.  But  with  the  LEEK  code  to  compute 
all  of  the  elastic  zones,  it  becomes  economically  feasible  to 
compute  the  problem. 

LEEK  is  a  very  sophisticated  elastic  code  and  may  be  used 
by  itself  as  a  powerful  tool  for  studying  seismic  wave  propaga¬ 
tion.  In  developing  it  we  have  made  novel  use  of  a  combination 
of  finite-dif fer  nee  techniques  commonly  used  by  seismologists 
and  by  physicists.  The  code  is  flexible  enough  to  be  able  to 
compute  problems  in  which  the  elastic  constants  are  arbitrary 
functions  of  position  and  in  which  zones  are  arbitrary  quadri¬ 
laterals.  In  addition,  either  stresses  or  displacements  may  be 
specified  as  the  boundary  conditions,  and  the  code  may  couple 
with  our  Lagrange  code  along  part  of  its  boundary.  A  brief 
description  follows  of  how  LEEK  works . 

Consider  a  linear  elastic  solid  in  cylindrical  coordinates 
about  the  z-axis.  Displacements  in  the  r  and  z  directions  art 
denoted  by  u  and  w,  and  because  of  the  (assumed)  cylindrical 

symmetry,  are  only  functions  of  r  and  z. 
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Axis  of 


Figure  A-l.  Zoning  for  crater  formation 
problem. 


The  equations  of  motion  are 
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A  finite-difference  mesh  of  quadrilateral  zones  is  defined 
which  covers  the  body.  Wherever  possible,  the  mesh  is  made  of 
uniform  rectangular  zones.  At  each  mesh  point  velocities  and 
displacements  are  defined.  Rectangular  "blocks"  of  mesh  points 
with  uniform  rectangular  zones  that  all  have  the  same  density 
and  elastic  constants,  X  and  U,  are  singled  out  for  a  fast 
finite-difference  acceleration  computation. 

Two  different  methods  are  used  to  compute  the  accelerations 
of  points.  Method  II  is  used  for  the  homogeneous  blocks  and  is 
described  below.  Method  I  is  used  for  all  other  points,  includ¬ 
ing  boundary  point3  and  interface  points  between  blocks,  except 
for  boundary  points  which  lie  '/ong  the  Lagrange-LEEK  interface. 
The  accelerations  of  these  points  are  computed  by  the  Lagrange 
code  from  forces  supplied  to  it  by  the  LEEK  code. 

In  Method  I,  direct  finite-difference  approximations  to 
Equations  (1)  through  (11)  are  used  to  compute  accelerations. 

The  finite-difference  equations  are  similar  to  those  described 
by  Wilkins  (Reference  A-l) .  In  this  scheme,  strains  are  first 
computed  in  all  existing  zones  surrounding  the  point  to  be 
accelerated.  Then  stresses  are  computed  from  Hooke's  law. 
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Stress  gradients  are  then  computed  at  the  grid  po-nt  ana  the 
point  is  accelerated  according  to  the  equations  of  motion.  This 

method  of  computing  accelerations  is  quite  general,  since  it 
allows  X,  g,  and  p  to  vary  from  zone  to  zone,  and  the  zones  need 

not  be  rectangular. 

A  much  faster  scheme,  Method  II,  is  used  to  compute  the 
accelerations  of  mesh  points  which  lie  within  homogeneous  rec¬ 
tangular  blocks.  In  this  case,  we  may  combine  Equations  (1) 
through  (11)  to  get  the  two  simpler  equations: 


where 


=  (X  +  2y)/p 

L 

4  =  p/p- 


The  subscripts  denote  partial  derivatives.  Simple  finite- 
difference  approximations  to  these  equations  can  be  made  (Refer 
ence  A-2) . 
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Once  the  accelerations  have  been  computed  by  either  Method  1 
or  Method  II.  they  can  be  integrated  over  a  t-e  step^to 
velocities  and  new  displacements.  Rlthoug  we 
omy  cylindrical  symmetry,  the  egun 
easily  derived  by  letting  r  -  ~  in  the  cynn 

equations . 
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APPENDIX  B 

DESCRIPTION  OF  THE  VARIOUS  MATERIAL  MODELS 


1.  YIELD  MODELS 

Two  types  of  yield  models  were  used  in  the  calculations — 
the  Mohr- Coulomb  and  the  von  Mises  models.  For  both  of  these 
models,  the  material  is  assumed  to  behave  elastically  up  to  some 
state  of  stress  at  which  yieldxng  occurs.  The  von  Mises  con¬ 
dition  for  flow  is  that  the  second  invariant*  of  the  stress 
deviator  tensor  is  a  certain  value  (Reference  B-l) .  For  an 
elastic-plastic  material,  the  same  result  obtains  when  volumetric 
compression  is  considered  to  be  purely  elastic,  and  when  there 
is  a  limit  to  the  amount  of  distortional  energy  that  can  be 
stored  elastically.  The  Mohr-Coulomb  condition  is  a  generali¬ 
zation  of  frictional  sliding.  The  shear  stress  required  for 
simple  slip  depends  on  the  cohesion,  c,  and  linearly  upon  the 
normal  pressure.  The  proper  generalization  for  this  condition 
under  any  state  of  stress  is  that  at  slip 

Y  =  c  +  bP 


*The  stress  de viator  tensor  is  the  stress  tensor  with  the 
average  stress  (pressure)  subtracted  from  the  diagonal  elements. 
The  second  invariant  of  a  tensor  of  rank  two,  t^j,  is  given  by 
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In  simple  tension, 


the  yield  stress  Y  -  >[3J2 


(Reference  B-2)  . 
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where  P  is  the  hydrostatic  component  of  stress  (Reference  B-3) . 

A  von  Mises  limit  is  conventionally  placed  on  the  yield  stress 
to  represent  ductile  flow  in  the  crystal  structure  as  an  upper 
limit  to  the  shear  strength.  For  both  of  these  yield  criteria, 
the  pressure  (and  volumetric  strain)  was  not  affected  by  plastic 

flow. 


2.  EQUATION  OF  STATE 

Two  models  of  hydrostatic  behavior  were  used  in  the  calcu¬ 
lations.  The  formula 


P  =  By 

represent  purely  elastic  behavior.  The  variable  y  is  a  com¬ 
pression  variable  useful  for  large  as  well  as  small  strains .  It 
is  computed  as 


y  = 


1, 


the  scaled  ratio  of  current  to  reference  densities.  The  co¬ 
efficient  B  is  the  (elastic)  bulk  modulus. 


The  other  equation  of  state  form  used  has  different  formulas 
for  loading  and  unloading 


p  -  v 


Loading 


P 


P  +  B„y  Unloading 
o  U 
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The  loading  modulus,  BL»  is  smaller  than  the  unloading  modulus, 
By.  The  figure  below  illustrates  the  calculational  procedure. 


Loading  to  high  pressure  (D  on  the  figure)  follows  the  path 

OECD,  where  OEC  has  the  slope  BL,  and  CD  has  the  slope  By.  Un 
loading  from  any  pressure  higher  than  C  to  zero  pressure  is  on 

the  line  DC A,  with  slope  By.  Any  further  loading  and  unloading 
is  on  the  line  PqACD.  If  the  loading  is  to  a  lower  pressure 
than  C,  the  crushing  is  not  complete.  One  such  loading  path  is 
OE.  Unloading  takes  place  along  EF  with  slope  By.  Reloading 

to  higher  pressure  from  F  takes  place  along  FECD. 
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